^A023479 


RADC-TR-76-55 

Final  Technical  Report 

March  1976 


COAT  MEASUREMENTS  AND  ANALYSIS 

Hughes  Research  Laboratories 


? 


Sponsored  by 

Defense  Advanced  Research  Projects  Agency 
ARPA  Order  1279 


Approved  for  public  release; 
distribution  unlimited. 


The  views  and  conclusions  contained  in  this  document  are  those  of 
authors  and  should  not  be  interpreted  as  necessarily  representing 
official  policies,  either  expressed  or  implied,  of  the  Defense 
Advanced  Research  Projects  Agency  or  the  U.  S.  Government. 


the 

the 


Rome  Air  Development  Center 
Air  Force  Systems  Command 
Grifflss  Air  Force  Base,  New  York  13441 


z' 


D D C 

SEmrisi 

APR  87  I97R 


in  1 1 « 1 1 III  TTci  U 


I 


This  report  has  been  reviewed  by  the  EABC  Infor^tlon  Offiee  (01)  and 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS 
it  will  be  releasable  to  the  general  public  including  fore  gn  na  lo 

This  report  has  been  reviewed  and  is  approved  for  publication. 


approved:  (Dr  JlJ 


ROBERT  F.  OGRODNIK 
Project  Engineer 


" iiiL  ' 'mt/in 


Do  not  return  this  copy. 


Retain  or  destroy. 


COAT  MEASUREMENTS  AND  ANALYSIS 

J.  E.  Pearson 
W.  P.  Brovm,  Jr. 

A.  Kokorowski 
M.  E.  Pedinoff 
C.  Yeh 


Co-ttractor : Hughes  Aircraft  Company 

Contract  Number:  F30602-75-C-0001 

Effective  Date  of  Contract:  2 July  1974 

Contract  Expiration  Date:  30  June  1975 

Amount  of  Contract:  $316,766.00 

Program  Code  Number:  5E20 

Period  of  work  covered:  Jul  74  - Jul  75 


Principal  Investigator: 
Phone : 


Mr.  J.  E.  Pearson 
213  456-6411, x283 


Project  Engineer: 

Phone : 


Robert  F.  Ogrodnik 
315  330-4306 


Approved  for  public  release; 
distribution  unlimited. 


This  research  was  supported  by  the  Defense  Advanced 
Research  Projects  Agency  of  the  Department  of 
Defense  and  was  monitored  by  Robert  F.  Ogrodnik 
(OCTM),  Griff iss  AFB  NY  13441. 


UNCLASSIFIED 


ICATIDN  DF  THIS  PAGE  (When  Dele  Entered) 


EPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2 GOVT  ACCESSION  NO 


3.  RECIPIENT'S  CATALOG  NUMBER 


y 


4.  TITLE  fand  Sufc(/(/eJ 


C^"L Measurements  and  Analysis# 


5.  TvpF  Of;  RFPnBT  A PFPinn  rnwFi|Fn 


Final  >Technical  /epi?t* 
1 Jui  74  - 1 Jul  75j 


'.E./f’earsorp  X'  M.E./Pedinoff  I 
W.P./trown,  Jr^  C./Yeh  A 


N/A 


^okorowski 


RFORUING  OT^GANIZATION  NAME  ANO  AOORESS 

Hughes  Research  Laboratories 
3011  Malibu  Canyon  Road 
Malibu.  CA  90265 


,8.  CONTRACT  OR  GRANT  NUMBERfs) 


F3/6;^2-75-C-40'^ 


■dm. 


^ JCCT.'T  ASk 

AREA  & WORK  UNIT  NUMBERS 


J.  0r;T'2790017 


II.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

Defense  Advanced  Research  Projects  Agency 
1400  Wilson  Blvd 
Arlington.  VA  22209 


no 


U.  MONITORING  AGENCY  NAME  8 ADDRESSfU  dllferml  from  Conirollint  DfflceJ 

Rome  Air  Development  Center  (OCTM) 

Griffiss  AFB.  NY  13441 


IS  SECURITY  CLASS,  (at  thU  reporlj 

UNCLASSIFIED 


ISe.  DECLASSIFICATION/ downgrading 
SCHEDULE 


JZL 


IS.  DISTRIBUTION  STATEMENT  fo/ 


Approved  for  public  release.  Distribution  Unlimited. 


17.  DISTRIBUTION  STATEMENT  (o!  l/ie  abalraci  aniaracf  in  Block  20.  U diUttmt  from  Report; 

Same 


18.  SUPPLEMENTARY  NOTES 

Robert  F.  Ogrodnik  Copies  are  available  in  DDC. 

OCTM 

AC315  330-4306 


19.  KEY  WORDS  (ConUnue  on  reverae  aide  If  neceaamry  end  Identify  by  bfock  nun  ber) 


CQATR'  Optical  Radar.  Laser  Phased  Arrays.  Thermal  Blooming  Compensation. 
Turbulence  Compensation.  Active  Optics.  Adaptive  Optical  Systems,  complex  target 
effects,  speckle  modulations 


20.  k^TRACT  (Continue  on  revarre  aide  If  neceeeery  end  identify  by  bfock  number) 


Coherent  optical  adaptive^^ techniques  (COAlX.have  been  studied  by  experiment, 
by  analysis,  and  by  computer  simulation.  The  experiments  have  utilized  a 21- 
channel.  visible-wavelength,  multidither  COAT  system,  while  the  computer  simula- 
tions have  dealt  with  both  multidither  (outgoing-wave)  and  phase-conjugate 
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i\20.  Abstract  (continued) 

the  status  of  the  21-channel  DARPA/RADC,  visible-wavelength,  multidither  COAT 
experimental  model  and  associated  hardware  and  its  use  in  the  experimental 
measurements  on  this  contract^ ^J^nropij^er  simulation  studies  of  propagation 
distortion  compensation  and  complex  tat*get  (speckle)  effects  are  also  described. 
Recoi^enilflj4.aQ&-Ara,made  for  future  work^ 

-MExperi mental  observations  with  the  21-channel  COAT  system  show  that  blooming 
distortions  occurring  in  the  first  30%  of  the  focused  propagation  path  can  be 
compensated,  leading  to  roughly  a factor  of  1.5  increase  in  peak  focused-beam 
i,jxadlj5in£ejw>  Depending  on  the  experimental  geometry,  however,  correction  factors 
l^from  1.0  to  4.0  have  been  observed.  Both  experimental  and  computer  simulation 
results  have  shown  that  turbulence  compensation  performance  is  not  degraded  by 
the  presence  of  thermal  blooming,  even  though  little  blooming  compensation  may 
occur.  Experimentally,  the  addition  of  tracking  and  focus  controls  to  18-channe 
planar-array,  COAT  phase  control  had  no  effect  for  blooming  distortions,  but 
produced  some  correction  for  an  artificially-generated  turbulence  alone.  The 
18-channel  phase  controls  could  remove  essentially  100, o of  the  turbulence  dis- 
tortions in  most  cases  where  there  was  no  significant  wavefront  tilt  error. 


omputer  simulation  of  phase-coi.jugate  (return-wave)  COAT  systems  has  shown 
almost  no  blooming  compensation.  \ In  addition,  there  is  some  evidence  that  a 
return-wave  COAT  system  may  redujje  the  focal -plane  irradiance  if  the  blooming 
distortions  become  too  strojis 

^Experimental  measurements  of  COAT  operation  with  equal  transmitter  and 
receiver  apertures  and  semidiffuse,  extended-glint  scotchlight  surfaces  have 
failed  to  produce  any  degradation  in  the  system  convergence  level,  even  though 
significant  spurious  modulations  within  the  dither  band  were  observed  in  the 
COAT  recei ver . B^xperimental  and  computer  simulation  results  have  shown,  however 
that  it  is  possible  to  reduce  the  COAT  system  convergence  level  if  spurious, 
multiplicative  r^eiver  noise  signals  are  strong  enough.  Such  signals  can  be 
produced  in  a COAP^receiver  when  the  speckle  pattern  return  (the  target  signatur^) 
moves  by  the  COAT  receiver.  A simple  analytical  model  has  been  developed  that 
can  accurately  predict  the  convergence  level  of  a multidither  COAT  system  in  the 
presence  of  spurious  receiver  modulation  noise.  These  studies,  as  well  as  other^, 
indicate  that  there  are  potential  problems  with  speckle-induced  noise  in  a COAT 
receiver.  Additional  studies  will  be  required,  however,  to  show  what  classes 
of  targets  and  scenarios  can  produce  spurious  receiver  signals  that  can 
significantly  degrade  the  performance  of  a multi  dither  COAT  system. 
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PREFACE 


This  report  was  prepared  by  Hughes  Research  Laboratories, 

Malibu,  California  under  Contract  F30602-75-C-0001 . It  describes  work 
performed  from  1 April  1975  to  1 June  1975  and  summarizes  the  principal 
results  of  the  entire  contract  during  the  period  1 July  1974  to  1 July  1975. 
The  principal  investigator  and  principal  scientist  is  Dr.  James  E.  Pearson. 
The  project  is  part  of  the  adaptive  optics  program  in  the  Opto- Electronics 
Department,  managed  by  Dr.  Viktor  Evtuhov,  at  the  Hughes  Research 
Laboratories. 
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glossary  OFAliURKVIATIONS  AND  ACRONYMS 


1 . 

COAT: 

Coherent  Optical  Adaptive  Techniques.  Any  one  of  a 
number  of  adaptive  optical  techniques.  Also  "sed  a^  a 
descriptive  modifier  as  in  "the  COAT  system." 

z. 

BAR: 

Beam  Active  Track.  Name  given  to  the  multidither 
cZtTci  controls  built  and  studied  on  this 

3. 

HRL: 

Hughes  Research  Laboratories,  Malibu,  California 

4. 

GRC: 

General  Research  Corporation,  McLean,  Virginia 

5. 

RADC: 

Rome  Air  Development  Center,  Rome,  New  York 

6 , 
7. 

NSWC; 

DARPA: 

Naval  Surface  Weapons  Center,  Silver  Spring,  Maryland 

Defense  Advanced  Research  Projects  Agency 
Washington,  D.C.  fe^^ncy. 

SUMMARY 


This  .-esearch  pro„ra™  o„  coherent  opt.cal  adaptive  techn.que. 

,COAT)  hae  three  primary  objectives:  (1)  to  experimentally  study  the 

u-d-ther  COAT  system  to  compensate  lor  thermal  coming 

ta-  m to  use  computer  simulations  to  cross-check  the  experi- 

memarrrsul'ts  and  to  investigate  areas  beyond  the  capabilities  of  the 

, , VO-  to  nerform  preliminary  studies  of  the  eueci 

experimental  hardware,  ( ) P TM=»r  for  mane  e. 

^ • r;  i=,tJon  on  COAT  system  convergence  and  perf 

of  speckle  noise  modulation  on  cu  y 

^ COAT  system  with  four  different  absorption  cells, 

uere  made 

two  liquid  and  two  gas.  ,,,  „„,,es 

structed  during  this  con  rac  . slightly  different  beam  geometry 

IRfiD  program.  Each  cell  was  used  in  si  ghtly 

as  discussed  in  Section  III,  but  each  geometry  can 

:„g  10.  b pm  scenario.  The  ratio  of  the  peak  V^T^or  I.  0 depend- 

—eli::";:  -:eoeeo  w!:e  0b:erved  between 

ino  on  the  experimental  geometry. 

ing  on  tne  greatest  correction 

blooming  obtained  with  liquid  or  gaseous  e t an??  of  the 

U observed  when  the  blooming  distortions  are  in  the  first  30%  of  the 

4-u  ikoaql-  correction  occurs  when  the  blooming 
focused  propagation  path;  the  rOAT  system 

r ad  to  the  last  50  to  60%  of  path.  In  all  cases,  the  COAT  system 

::  enc  llvel  was  stable,  with  no  evidence  of  oscillation  or  instabil  ty 

Ttbe  Loming  distortions  were  increased  well  beyond  that  present  with 

the  optimum  transmitter  power. 


PRECEDI  BU-NK-NOT  flM'oD 


15 


Multidither  adaptive  tracking  and  foe  (S  controls  were  designed  and 
constructed  to  be  used  with  the  18-channel  COAT  system.  The  design 
goals  of  a 100  Hz  tracking  sy  ..em  and  a iO  Hz  focusing  system  were  m.et  by 
the  hardware.  In  the  thermal  blooming  measurements,  these  controls  pro- 
duced no  further  increase  in  peak  target  irradiance  over  that  provided  by 
the  18-channel  plana.r  array;  in  one  turbulenc('-only  experiment,  however, 
some  correction  was  observed. 

We  have  experimentally  studied  the  effects  of  target  motion- induced 
speckle  modulation  on  the  18-channel  COAT  system.  A scaled  laboratory 
experiment  using  a rotating  scotchlite  target  produced  a wideband  speckle 
noise  spectrum  over  the  dither  frequency  band  which  did  not  observably 
interfere  with  the  convergence  properties  of  the  system.  In  another  series 
of  experiments,  an  acousto-optic  modulator  with  64%  modulation  efficiency 
was  placed  in  the  DARPA/RADC  COAT  optical  system  and  swept  over  a 
variable  frequency  bandwidth  in  a time  comparable  to  the  servo  loop  con- 
vergence time.  The  speckle  modulation  artificially  generated  by  this  tech- 
nique produced  interference  signals  over  the  dither  band  greater  than  the 
normal  dither  signals,  and  produced  converged  power  degradation  of  35%. 

Computer  simulation  modeling  of  COAT  system  performance  in  the 
presence  of  speckle  noise  modulation  has  been  used  to  obtain  the  average 
convergence  level  of  an  18-channel  multidither  COAT  system  in  the  presence 
of  this  noise.  An  oscillation  of  the  system  about  the  average  convergence 
level  has  been  observed.  A statistical  analysis  of  the  COAT  system  and  of 
the  speckle  noise  has  been  developed  for  use  in  explaining  the  simulation 
results.  Good  agreement  has  been  observed  between  the  analytical  model 
and  the  computer  simulation  data.  The  experimental  and  theoretical 
results  obtained  thus  far  indicate  that  the  system  performance  can  be 
significantly  degraded  if  the  modulation  is  large  enough,  but  the  systematic 
evaluation  of  which  targets,  target  ranges,  and  target  motions,  if  any,  can 
actually  produce  such  large  modulations  are  being  pursued  on  other  pro- 
grams. To  date,  however,  our  experimental  results  with  the  18-channel 
COAT  system  (equal  transmitter  and  receiver  apertures)  and  realistic 
targets  indicate  that  speckle  modulations  will  not  seriously  degrade  COAT 
system  performance.  This  is  a significant  conclusion  since  any  operational 
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CO/iT  system  will  deal  with  semidiffuse,  multiple  glint  targets  that  do 
produce  speckle-like  target  signatures. 

Turbulence  compensation  in  the  presence  of  thermal  blooming  has 
been  studied  by  experiment  and  by  computer  simulation.  We  have  observed 
that  turbulence  compensation  is  not  significantly  degraded  by  the  existence 
of  thermal  blooming,  even  though  the  thermal  blooming  distortions  may  not 
be  removed  by  the  COAT  system.  This  result  is  particularly  significant  for 
laser  systems  operating  at  3.  8 pm  or  at  10.  6 pm  with  larger  (>1  m)  aper- 
tures where  turbulence  distortions  dominate  the  beam  degradation.  Just  as 
blooming  and  turbulence  can  be  treated  separately  and  later  combined  to 
determine  the  total  beam  distortion,  so  can  the  COAT  compensation  for  each 
type  of  distortion  be  treated  separately.  In  the  experimental  measurements, 
identical  compensation  performance  was  observed  with  single  glint  targets 
and  with  semidiffuse,  rough-surface,  rotating  spherical  targets. 

Our  computer  simulation  studies  of  COAT  compensation  for  thermal 
blooming  have  been  restricted  to  return-wave  (e.  g.  , phase-conjugate)  COAT 
systems  because  of  the  large  computational  cost  when  the  multidither  servo 
code  is  combined  with  the  nonlinear  propagation  code.  This  point  is  dis- 
cussed further  in  Section  IV.  The  operation  of  return-wave  systems  is  sig- 
nificantly different  with  nonli.near  distortions  than  that  of  a multidither, 
outgoing- wave  system  (like  the  18-channel  experimental  system).  We  have 
observed  very  little  blooming  compensation  in  the  simulations  and  have 
found  that  return-wave  systems  can  actually  reduce  the  target  irradiance 
when  the  transmitted  power  exceeds  the  optimum  power  (that  power  which 
produces  maximum  target  irradiance).  We  suspect  that  this  type  of 
behavior  may  be  characteristic  of  any  return-wave  system,  including  those 
that  employ  compensated  imaging  concepts  (imaging-COAT  systems). 

Based  on  previous  data  ''  and  on  the  investigations  summarized  in 
this  report,  we  conclude  that  a multidither,  outgoing- wave  COAT  system 
can  remove  70  to  100%  of  turbulence -induced  beam  distortions,  100%  of 
fixed  beam  distortions  (optics,  etc.),  can  provide  about  a factor  of  1 . 5 to 
2.0  improvement  in  peak  beam  irradiance  in  the  presence  of  thermal 


■•'J.  E.  Pearson,  "Coherent  Optic."xl  Adaptive  Techniques,  " Contract  F30602- 
73-C-0248m  Final  Technical  Report,  RADC-TR-  75-46.,  Jan.  1975. 
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blooming,  and  can  accomplish  all  of  these  functions  simultaneously,  limited 
only  by  the  number  of  correction  channels  and  by  the  system  speed  of 
response.  The  compensation  can  be  achieved  with  realistic  targets  having 
multiple  glints,  or  semidiffuse  surfaces.  Propagation  of  high  energy  laser 
beams  through  the  atmosphere  to  produce  large  focal-plane  irradiances 
should  thus  not  be  limited  by  the  atmosphere  if  suitable  COAT  system 
designs  are  used  and  if  severe  thermal  blooming  distortions  can  be  avoided. 
Additional  work  is  indicated  in  several  areas:  (1)  speckle  modulation  effects 
to  determine  what  targets  and  scenarios  can  produce  speckle  modulations 
that  are  strong  enough  to  degrade  multidither  COAT  performance  and  to 
design  COAT  systems  that  are  insensitive  to  such  modulations;  (2)  thermal 
blooming  compensation  — to  investigate  techniques  that  can  minimize  thermal 
blooming  distortions;  (3)  computer  simulation  — to  develop  codes  that  can 
simulate  thermal  blooming  compensation  with  multidither  COAT  systems 
and  realistic,  multiple- glint  targets;  (4)  COAT  operation  at  infrared  wave- 
lengths — to  demonstrate  the  hardware  technology  required  for  COAT  sys- 
tems at  longer  wavelengths  and  higher  powers  and  to  investigate  target 
effects  (speckle,  glint  destrxiction)  with  realistic  target  surfaces  and 
geometries;  (5)  new  COAT  concepts  — to  investigate  experimentally  new 
servo  and  system  concepts  that  offer  improved  compensation  performance, 
greater  versatility,  greater  simplicity  in  hardware  design  and  implementa- 
tion and  more  reliable  operation  under  a wider  variety  of  operational 
conditions . 
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INTRODUCTION 


A.  “Program  Objectives 

There  are  two  primary  objectives  of  this  program.  The  first 
objective  is  to  determine  the  performance  limits  of  multidither  coherent 
optical  adaptive  techniques  (COAT)  by  performing  scaled  laboratory  experi- 
ments which  are  designed  to  produce  quantitative  data  on  the  nature  of 
thermal  blooming  and  turbulence  and  on  the  ability  of  multidither  COAT  to 
correct  for  these  distortions,  particularly  when  the  target  has  many  moving, 
time-varying  glints.  The  second  objective  is  to  use  computer  simulation  as 
a cross-check  on  the  experiments  and  as  an  analytical  tool  to  extend  the 
understanding  of  propagation  distortions  and  of  COAT  systems  beyond  the 
capabilities  of  the  experimental  hardware. 

B.  Research  Program  Plan 

The  performance  of  multidither  COAT  has  already  been  proved  with 
an  18- element  visible  wavelength  system  developed  and  tested  on  DARPA/ 
RADC  Contract  F30602-73-C-0248  which  was  concluded  in  July  1974.  This 
same  system  is  used  for  the  laboratory  experiments  in  this  program.  The 
computer  simulation  codes  for  atmospheric  turbulence,  thermal  blooming, 
and  the  COAT  system  were  developed  on  DARPA/NOL  Contract  N60921- 
74-C-0249.  The  previous  DARPA/RADC  COAT  contract  also  supplied  a 
design  for  a gas  absorption  cell  which  has  been  used  to  simulate  convection- 
dominated  thermal  blooming  in  the  atmosphere. 

The  research  program  for  this  contract,  illustrated  in  Fig.  1,  runs  . 
from  1 July  1974  through  30  June  1975.  The  required  oral  presentation  was 
made  on  3 December  1974  as  part  of  the  DARPA/NOL  Adaptive  Optics 
Symposium  held  at  Lincoln  Laboratories.  A contract  amendment  negotiated 
during  the  second  contract  quarter  provided  for  the  addition  of  auto-tracking 
and  autofocus  controls  to  both  the  DARPA/RADC  COAT  hardware  and  the 
Naval  Ordnance  Laboratories  computer  code.  A second  amendment,  nego- 
tiated during  the  third  quarter,  eliminated  a high-power  design  guidelines 
task  so  that  time  and  funds  could  be  directed  toward  experiments  on  COAT/ 
target- signature  interactive  effects.  The  schedule  in  Fig.  1 reflects  these 
changes. 
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This  report  serves  two  functions:  (1)  it  presents  the  accomplishments 
of  the  fourth  contract  quarter  (1  April  1975  to  1 July  1975);  (2)  it  summarizes 
the  principal  accomplishments  of  the  entire  contract.  The  report  is  thus  both 
the  fourth  quarterly  technical  report  and  the  final  contract  report. 

Section  II  reviews  the  status  of  the  experimental  hardware  thkt  has 
come  to  be  known  as  the  DARPA/RADC  COAT  system.  Items  summarized 
include  the  18- channel  planar  COAT  array,  the  3-channel  Beam  Active 
Track  (BAT)  system,  the  associated  optical  systems  used  in  experimental 
laboratory  studies,  and  the  laboratory  distortion  generators  used  to  produce 
thermal  blooming  and  turbulence. 

Section  III  summarizes  the  experimental  results  obtained  on  turbu- 
lence and  thermal  blooming  compensation.  Preliminary  investigations  on 
COAT  system  operation  with  complex,  extended-glint  targets  are  presented. 
The  section  briefly  discusses  the  effects  of  spurious  receiver  modulations 
produced  by  moving  speckle-patterns  (the  target  signature). 

Section  IV  summarizes  the  results  of  the  analytical  efforts  directed 
toward  investigating  turbulence  and  thermal  blooming  compensation  as  well 
as  complex  target  effects.  The  data  presented  were  obtained  primarily  by 
computer  simulation  of  the  propagation  distortions  and  the  COAT  system 
(either  phase  conjugate  or  multidither). 

Section  V summarizes  the  principal  results  of  the  contract  and  the 
conclusions  that  can  be  drawn  from  this  work.  The  implications  for  addi- 
tional research  that  should  follow  this  work  are  presented. 

r ^ 
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II. 


COAT  SYSTEM  HARDWARE 


A.  Planar  Array  Phase  Controls 

The  original  DARPA/RADC  COAT  system  was  developed  on 

Contract  F30602-73-C-0248.  A detailed  description  of  the  system  is  pre- 

1-5 

sented  in  the  five  contract  reports  " for  that  program.  A block  diagram 
of  the  complete  system  is  shown  in  Fig.  2 and  Fig.  3 is  a photograph  of  the 
signal  conditioning  and  servo  electronics  for  the  18-phase  control  channels. 
The  features  and  performance  of  the  system  are  summarized  in  Tables  1 
and  2.  As  can  be  seen  from  Table  2,  all  of  the  fundamental  multidither 
COAT  concepts  have  been  successfully  demonstrated  with  this  system.  The 
remaining  issues  are  no  longer  concerned  with  whether  or  not  the  system 
will  work,  but  rather  with  how  well  it  will  work  under  various  conditions. 
The  data  obtained  on  this  contract  answer  the  latter  question  for  some 
cases  and  succeeding  programs  now  in  progress  are  designed  to  further 
quantify  the  degree  of  successful  COAT  compensation  for  turbulence  and 
blooming  with  realistic  targets. 


Table  1.  DARPA/RADC  Planar  Array  COAT  System  Features 
1.8  W,  0.488  pm,  argon  laser 

Planar  output  arrays  with  any  configuration  containing  up  to  18  phase- 
controlled  elements  and  one  reference  element. 

PZT  cylinders  and  bimorphs  used  as  phase  dither  and  phase  correc- 
tion drivers  ("piston"  elements).  Up  to  60°  peak  dither  modulation 
and  over  ±2  \ correction  available  at  X.  = 0.  488  pm. 

Dither  frequencies  adjustable  from  8 to  32  kHz;  1 channel/frequency 
or  sine/cosine  operation. 

Chopper- stabilized  AGC  with  46  dB  dynamic  range  and  1 msec 
response  time. 

Sample-and-hold,  slew-offset  capability  — electronic  and  mechanical. 

Two  photomultiplier  detectors  available;  one  for  "local  loop"  opera- 
tion (no  path  distortion),  one  for  "target  loop"  operation  (path  contains 
turbulence  and  blooming). 
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COAT  system  block  diagram  showing  additional  multidither  servo 
loops  for  tracking  and  focus  control. 


RADC/COAT  multidither  electronics:  18-phase- 
control  servo  channels  are  shown;  the  2 tracking 
(x  and  y)  channels  and  1 focus  channel  are  shown  in 


Table  2.  DARPi\/RADC  Planar  Array  Performance 


Convergence  time  of  1 . 5 ms  or  less  (10%  to  90%). 

Good  agreement  between  computer  simulation  of  servo  and  experi- 
mental results. 

Single  glint  tracking  at  rates  up  to  30  mrad/ sec. 

1.3  X diffraction-limited  performance  (no  propagation  distortions). 
Convergence  on  strongest  glint  in  a moving,  multiple  glint  target. 
Demonstrated  turbulence  compensation. 

Demonstrated  thermal  blooming  compensation. 

Demonstrated  turbulence  compensation  in  the  presence  of  blooming. 

Demonstrated  compensation  performance  with  "extended-glint" 
targets. 


B.  Multidither  Tracking  and  Focus  Controls 

Amendment  No.  1 to  the  original  contract  provided  supplemental 
funds  to  add  three  additional  multidither  control  channels  to  the  18 -channel 
planar  array  controls:  2 tracking  controls  and  1 spherical-focus  control. 

Although  we  usually  will  refer  to  the  entire  COAT  system  as  the  "21 -channel 
DARPA/RADC  COAT  system,"  we  will  have  occasion  to  refer  to  just  the 
focus  and/or  tracking  controls.  In  this  report,  we  will  use  the  term 
"18-channel"  to  refer  to  the  18-element  planar  array  system  which  pro- 
duces piston-type  phase  control  on  18  separate  transmitted  beam  segments. 
The  acronym  "BAT"  (for  Beam  Active  Tracking)  will  be  applied  to  the 
focus  and/or  tracking  controls  alone.  When  sp  .aking  of  combinations,  we 
will  use  the  terms  "20-channel"  (COAT  plus  tracking),  or  "21-channel" 
(COAT  plus  focus  plus  tracking). 

The  tracking  and  focus  controls  are  shown  in  Fig.  4 and  Fig.  5 
is  a schematic  of  the  focus  actuator.  The  servo  control  channels  were 
built  by  dividing  the  servo  parameters  of  the  18-channel  system  by  10  for 
the  tracking  channels  and  by  100  for  the  focus  control.  This  procedure 
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was  expected  to  produce  a tracking  servo  with  a 1 kHz  dither  and  a 10  msec 
convergence  time  (100  Hz  response)  and  a focus  servo  with  a 100  Hz  dither 
and  a 100  msec  convergence  time  (10  Hz  resoonse).  The  oscilloscope  traces 
of  target  irradiance  versus  time  shown  in  Fig.  6 demonstrate  the  achieve- 
ment of  these  design  goals. 

Both  controls  could  be  made  to  operate  faster  since  the  lowest 
resonance  frequency  of  the  focus  actuator  is  3.4  kHz  and  the  tracking 
galvanometer-plus-mirror  resonance  is  at  1.7  kHz.  Although  we  chose  to 
operate  the  two  tracking  controls  with  two  separate  dither  frequencies 
(750  Hz  and  1 kHz),  sine-cosine  operation  will  work  equally  well  (producing 
a conical-scan  tracker). 

The  tracking  performance  of  the  controls  is  shown  in  Fig.  7.  An 

interesting  behavior  is  observed:  the  controls  do  not  automatically  steer 

the  beam  boresight  axis  onto  the  glint  when  the  phase  and  BAT  controls  are 

used  together.  The  stable  lock  position  depends  on  the  initial  conditions 

as  shown  in  the  figures.  This  behavior  is  not  entirely  reproducible  and  is 

a consequence  of  the  sidelobe  structure  of  the  segmented  transmitting 

aperture  and  of  the  choice  of  dither  frequencies.  (The  behavior  would  not 

exist  if  the  tracking  servo  were  faster  than  the  phase  control  servo.)  This 

explanation  is  discussed  in  greater  detail  in  the  third  contract  quarterly 
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report . 

Although  the  focus  control  servo  has  performed  as  predicted  by  the 
design,  the  actuator  shown  in  Fig.  5 has  one  problem:  it  in:  -oduces  a 

small,  but  not  insignificant  amount  of  aberration  into  the  optical  beam. 

When  the  actuator  is  used  with  the  IS-channel  COAT  system,  the  distortion 
is  no  problem  since  the  COAT  system  can  completely  remove  it.  The  dis- 
tortion does  not  allow  the  actuator  to  be  used  by  itself,  however,  without 
introducing  undesired  beam  distortions.  Consequently,  all  of  the  experi- 
ments discussed  in  Section  III  that  involve  the  focus  control  also  use  the 
1 8 - channel  COA  T system. 

The  BAT  tracking  and  focus  controls  have  also  been  implemented 
into  a multidither  computer  simulation  that  uses  a deformable  mirror  as  a 
dither  and  corrector  element.  The  nonboresighting  effects  seen  in  tlie 
experiment  were  thus  not  observed.  We  also  decided  to  implement  a 
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Fig.  7 . 

Tracking  performance  of  COAT  controls,  (a)  Tracking  control 
off.  The  formed  array  pattern  is  also  shown  for  reference. 

(b)  Tracking  control  performance.  Curves  with  tracking  control 
only  and  with  both  phase  and  tracking  controls  are  shown.  The 
ability  of  the  system  to  lock  in  a nonoptimum,  off-boresight  con- 
dition is  evident. 


two-axis  focus  control  in  the  simulation,  in  effect  using  two  cylindrical  lens 
controls  that  focus  along  orthogonal  axes.  This  type  of  focusing  should  be 
particularly  useful  against  thermal  blooming  since  the  negative  thermal 
lens  is  much  stronger  perpendicular  to  the  wind  direction  than  parallel 
with  it. 

The  computer  simulation  results  have  shown  the  expected  10  msec 
convergence  time  for  the  tracking  controls  and  100  msec  for  the  focus. 

With  the  focus  control,  however,  we  observed  another  initially  unexpected 
result:  the  final  converged  level  can  depend  on  the  magnitude  of  the  initial 
defocus.  This  behavior  is  now  understood  to  be  a generalized  manifestation 
of  the  "ZnTT  problem"  encountered  with  deformable  mirror  COAT  systems. 
This  simulation  has  been  used  to  study  turbulence  and  blooming  compensa- 
tion; results  of  these  calculations  are  discussed  in  Section  IV, 

C . Thermal  Blooming  Generators 

Several  different  absorption  cells  have  been  built  during  the  course 
of  this  contract  and  each  one  has  been  used  in  a variety  of  experimental 
arrangements.  Both  liquid  and  gaseous  media  have  been  used.  We  did  not 
observe  any  medium-dependent  differences  in  the  nature  of  blooming  or 
the  amount  of  COAT  correction  that  could  be  produced. 

1 . Gas  Cells 

5 

Several  designs  were  evaluated  for  a cell  to  produce 
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convection-dominated  thermal  blooming.  The  first  design  we  wmployed  is 
shown  in  Fig.  8.  The  design  uses  a stationary  gas  cell  (Fig.  8(a))  and  a 
moving  mirror  arrangement  to  simulate  a transverse  wind  (Fig.  8(b)). 
Experiments  in  which  the  cell  was  mounted  vertically  to  minimize  convec- 
tion effects  gave  substantially  the  same  results  as  those  in  which  the  cell 
was  horizontal  (see  Section  III  of  this  report). 

5 

The  initial  cell  design  identified  SF^,  seeded  with  NO^  for  absorp- 
tion, as  the  best  gas  for  producing  maximum  blooming  effects.  For  reasons 
that  are  still  not  apparent,  SF^  did  not  produce  the  expected  amount  of 
blooming.  We  theorize  that  the  large  number  of  V-T  states  in  SF^  pre- 
vents rapid  transfer  of  the  energy  absorbed  by  the  NO^  into  thermal 
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Statis  gas  absorption  cell,  (a)  Photograph  of  cell. 

(b)  Schematic  of  "wind"  generation  apparatus.  The  moving 
mirror  causes  the  input  beam  to  move  relative  to  the  gas 
cell.  The  fixed  mirrors  bring  the  moving  beam  to  the 
back  of  the  moving  mirror  so  that  the  output  beam  is 
stationary. 


(translational)  energy  in  the  SF^.  The  next  choice  for  a buffer  gas  was  Xe 
and  it  produced  blooming  according  to  our  predictions.  The  disadvantage 
of  Xe  is  its  much  higher  cost  compared  with  SF^. 

Although  the  static  gas  cell  in  Fig.  8 operates  satisfactorily,  it  has 
two  serious  drawbacks.  First,  it  cannot  be  used  with  slewing  beams  because 
of  limited  input  and  output  aperture  sizes.  More  important,  however,  the 
optical  arrangement  places  the  blooming  distortions  in  the  last  55%  of  the 
focused  propagation  path.  This  is  not  only  an  unrealistic  propagation 
scenario,  but  as  discussed  in  Section  III,  it  precludes  observing  any  sig- 
nificant COAT  compensation  for  blooming.  To  correct  these  deficiencies, 
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we  built  a flowing  gas  cell  ' as  shown  in  Fig.  9. 

The  flowing  gas  cell  can  produce  transverse  wind  velocities  up  to 
12  cm/sec  as  measured  by  a hot-wire  anemometer.  The  fill  gases  are 
again  Xe  and  NO^.  No  turbulence  effects  are  observed  indicating  good 
laminar  flow.  This  cell  was  used  in  the  same  vertical  configuration  as 
the  static  gas  cell  to  minimize  buoyancy  effects  although  our  observations 
have  indicated  that  this  is  an  unnecessary  complication  of  the  optics.  The 
optical  arrangement  is  similar  to  Fig.  8(b)  except  that  the  focusing  lens 
can  be  placed  much  closer  to  the  cell  input  window  and  the  mirror  labeled 
1 or  2 in  Fig.  8(b)  is  stationary. 

We  experienced  some  difficulties  in  using  the  flowing  gas  cell 
mounted  as  shown  in  Fig.  9(b).  First,  the  blower  was  not  vacuum  tight 
so  that  small  quantities  of  NO^  slowly  leak  out  as  it  runs.  The  leak  is 
too  small  to  cause  a safety  hazard,  but  it  is  annoying  because  of  the 
unpleasant  smell  of  N02*  Second,  the  optics  were  cumbersome  and 
numerous  and  had  just  enough  residual  distortion  to  prevent  accurate 
measurement  of  the  beam  distortion  produced  by  blooming  alone  with  no 
COAT  correction.  Since  natural  convection  effects  are  not  noticeable 
for  our  experimental  conditions,  a horizontal  mounting  arrangement  could 
be  used.  This  change  will  probably  eliminate  the  optics  problem,  but  will 
not  help  the  NO^  leak  problem.  Time  on  this  contract  lid  not  permit  us  to 
try  out  a horizontal  mounting,  but  we  expect  to  use  such  an  arrangement 
in  future  planned  experiments. 
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(b)  Cell  is  mounted  with  a vertical  path  to 
minimize  buoyancy  effects. 


Flowing  gas  cell  for  blooming  studies 


2. 


Liquid  Cells 

All  of  the  thermal  blooming  studies  on  this  program  have 
been  scaled  to  simulate  conditions  encountered  by  focused  beams  propagating 
over  relatively  short,  low  altitude  paths  in  the  atmosphere.  The  scaling 
considerations,  discussed  in  detail  in  Section  III-A,  assume  that  the  heat 
transfer  is  dominated  by  forced  convection  (wind),  which  is  the  normal  con- 
dition in  the  real  atmosphere.  As  long  as  this  assumption  is  satisfied,  and 
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if  we  are  not  concerned  with  effects  like  convective  cooling  that  occur  only 
at  high  altitudes,  there  are  no  fundamental  differences  between  liquid  and 
gaseous  blooming  media.  The  differences  that  have  been  observed  occur 
only  for  collimated-beam  propagation  and  not  for  focused  beams  such  as  we 
are  considering. 

Liquid  media  have  two  distinct  advantages  over  gases  for  our  visible- 
wavelength  experiments.  First,  a flowing-liquid  cell  is  simple  to  construct 
and  the  liquids  are  nontoxic  and  easy  to  handle.  Second,  much  stronger 
blooming  effects  can  be  produced  for  a given  optical  power  and  absorption 
because  liquids  are  available  with  large  values  of  dn/dT,  the  temperature 
coefficient  of  the  refractive  index. 

As  part  of  the  Hughes  Research  Laboratories  IR&D  program,  we 
have  studied  COAT  compensation  for  thermal  blooming  in  liquids.  The 
results  obtained  are  included  with  this  contract  work  because  of  their 
applicability  to  the  questions  being  asked  on  this  program  and  because  the 
DARPA/RADC  COAT  system  was  used  in  the  measurements.  Two  liquid 
cells  have  been  used:  a 1 mm  thick  static  cell,  and  a 20  cm  long  flowing 
liquid  cell.  The  flowing  cell  is  shown  in  Fig.  10.  We  have  used  carbon 
tetrachloride  (CCl^)  methanol  for  liquids  with  iodine  added  to  absorb  the 
0.  488  pm  laser  light.  Methanol  is  more  convenient  to  use  than  CCl^  in  the 
flowing  cell  because  the  very  high  value  of  dn/dT  in  CCl^  makes  it  hard  to 
produce  a uniform,  turbulence-free  flow. 

D.  Turbulence  Phase  Screen 

Generating  artificial  turbulence  in  the  laboratory  is  not  difficult,  but 
doing  it  in  a calibrated  fashion  is.  We  chose  to  use  a single  phase  screen 
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(b) 

Fig.  10.  Flowing  liquid  cell  for  thermal  blooming  studies. 

(a)  Schematic  diagram,  (b)  Photograph  of  20  cm  cell. 
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which  is  placed  just  ahead  of  the  lens  which  focuses  through  the  blooming 
cell.  The  phase  screen  is  fabricated  as  shown  in  Fig.  11(a)  by  sputtering 
patches  of  SiO^  onto  a fused  quartz  plate  through  a mask.  The  mask  is  a 
thin  metal  plate  with  equal  size  holes  drilled  in  it  in  a regular  pattern. 

O O 

Three  to  six  evaporations  of  about  1200  A of  SiO^  (X/4  at  4880  A)  are  made 
with  the  mask  in  a different  position  each  time.  In  this  way  the  plate  is 
nearly  covered  with  a random  patch  pattern  and  the  phase  variations  range 
from  zero  to  as  much  as  3X./2.  An  electronic  micrometer  scan  through  one 
of  these  phase  plates  is  shown  in  Fig.  11(b).  The  peak-to-peak  phase 
variation  is  about  3/4  wavelength  and  the  minimum  size  of  a given  "turbule" 
is  about  0.  5 mm.  The  size  of  an  element  in  the  COAT  array  is  also  about 
0.  5 mm  at  the  location  where  the  plate  is  used.  To  produce  time-varying 
distortions,  the  plate  is  rotated  through  the  beam  as  shown  in  Fig.  11(b). 

This  type  of  turbulence  phase  screen  has  two  disadvantages.  First, 
the  turbulence  occurs  only  in  one  plane,  unlike  atmospheric  turbulence 
which  is  a distributed  path  effect.  Second,  the  magnitude  of  the  distortions 
and  their  spatial  frequency  spectrum  are  quite  different  from  atmospheric 
turbulence.  In  particular,  we  have  found  that  our  phase  screens  contain 
very  little  tilt  or  beam  steering  distortion;  tilv  is  one  of  the  larger  compo- 
nents of  atmospheric  turbulence.^  ^ 

The  artificial  nature  of  the  turbulence  is  not  important  for  the  pur- 
poses of  this  program,  however.  We  have  already  demonstrated  the  ability 

512  813 

of  COAT  to  correct  for  turbulence  ’ and  we  have  shown  analytically  ' 

that  turbulence  and  blooming  effects  can  be  considered  separately  and 
then  combined  later  to  give  an  effective  overall  reduction  in  beam  Strehl 
ratio.’''  In  this  contract,  we  are  investigating  the  correction  for  turbulence- 
like effects  in  the  presence  of  thermal  blooming. 

E..  Laboratory  Optical  Arrangement 

The  details  of  the  optical  layout  differ  for  each  experiment.  Each 
arrangement  has  several  common  features,  however,  and  in  functional 


■i'Strehl  ratio  is  the  ratio  of  the  peak  beam  irradiance  to  the  diffraction- 
limited  irradiance. 
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Fig.  11.  Phase  screen  for  generating  artificial  turbulence 

effects,  (a)  Evaporated  random  patches,  Si02  on  a 
fused  quartz  substrate,  (b)  Calibration  scan  through 
an  artificial  turbulence  plate  using  an  electronic 
micrometer.  A change  in  thickness  of  4880  K will 
produce  a one  wavelength  phase  shift  for  our  experi- 
ments. (c)  Time  variations  produced  by  rotating  plate 
through  the  laser  beam. 
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block  form  they  all  look  like  the  schematic  in  Fig.  12.  The  first  reducing 
telescope  recollimates  the  transmitting  array  to  a size  usable  by  the  BAT 
actuators.  There  are  two  optical  paths  that  have  identical  optics  and  differ 
only  in  that  one  path  has  the  turbulence  and  blooming  distortions  in  it.  The 
second  undistorted  path  is  required  so  that  the  COAT  system  can  be  used  to 
remove  random  step  phase  differences  introduced  by  the  "phasor  matrix." 

When  the  photomultiplier  (PMT)  at  location  No.  2 is  connected  to  the 
COAT  servo  system,  a diffraction-limited  beam  is  put  into  the  path  contain- 
ing the  BAT  actuators  and  the  distortions.  The  Strehl  ratio  with  no  COAT 
correction  is  measured  this  way.  The  degree  of  COAT  correction  is  mea- 
sured bv  switching  the  PMT  at  location  No.  1 into  the  servo  system.  This 
comparison  procedure  is  possible  only  if,  when  the  blooming  and  turbulence 
distortions  are  removed,  there  is  no  change  if  the  beam  at  the  primary 
targf  t plane  as  the  PMT  choice  is  switched  from  No.  2 to  No.  1 and  back. 

As  mentioned  previously  and  as  discussed  further  in  Section  III,  this 
requirement  has  limited  our  ability  to  measure  the  effectiveness  of  focus 
control  alone  in  compensating  turbulence  and  blooming. 

For  all  the  experiments  discussed  in  Section  III,  we  will  use  the 
following  terminology!  (1)  "NO  COAT  correction"  means  PMT  No.  2 in 
Fig.  12  is  used;  (2)  "COAT  Corrected"  means  PMT  No.  1 is  used.  Occa- 
sionally, we  may  refer  to  "local  loop  back,"  meaning  use  of  PMT  No.  2, 
or  "target  loop  lock,"  meaning  use  of  PMT  No.  1.  In  all  cases,  the  beam 
profile  and  the  Strehl  ratio  are  measured  at  the  "primary  target  plane" 
indicated  in  Fig.  12. 


".ajai.. 


III.  EXPERIMENTAL  TECHNICAL  ACCOMPLISHMENTS 

The  experimental  work  on  this  p'-ogram  consists  of  five  major  tasks; 
(1)  construction  of  a gas  cell  thermal  blooming  simulator,  (2)  construction 
of  multidither  tracking  and  focus  controls,  (3)  measurement  of  COAT  com- 
pensation for  thermal  blooming  with  single  and  multiple  glint  targets, 

(4)  measurement  of  COAT  compensation  for  turbulence  and  thermal  bloom- 
ing simultaneously  with  complex  targets,  (5)  preliminary  study  of  target- 
signature  / speckle  effects  on  a multidither  COAT  system.  The  last  three 

tasks  also  have  couterparts  in  the  analytical  studies  discussed  in  Section  V. 

^ 6 8 

The  previous  section,  as  well  as  earlier  reports,  summarized  the  per- 
formance of  the  hardware  that  was  built.  This  section  summarizes  the 
results  of  the  experimental  measurement  tasks, 

A . Thermal  Blooming  Compensation 

1 . Scaling  Considerations 

The  thermal  blooming  experiments  on  this  program  are 
designed  to  increase  our  understanding  of  thermal  blooming  and  of  the 
limitations  of  multidither  COAT  systems  in  trying  to  remove  these  beam 
distortions.  The  experiments  are  thus  scaled  so  that  the  results  can  be 
extended  to  other  wavelengths,  ranges,  etc.  Since  blooming  effects  are 
most  pronounced  at  10.  6 fim,  we  chose  to  scale  our  visible  wavelength 
laboratory  experiments  to  a typical  10.6  pm  scenario.  In  order  to  scale  a 
10.6  pm,  long-path  experiment  to  a visible  wavelength  laboratory  e..peri- 
ment,  three  things  must  be  considered:  (1)  wind  speed  and  beam  SiZe  to 

ensure  convection-dominated  heat  transfer  (the  only  case  of  interest  with 
actual  atmospheric  conditions),  (2)  total  absorption,  and  (3)  blooming 
strength  (transmitter  power). 

The  wind  speed  and  beam  size  are  related  to  the  medium  properties 
by  the  Peclet  number  defined  by 
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where  V is  the  transverse  winrl  velocity,  r^  is  ttie  l^e.im  i-ririius  ;ind  the 
medium  properties,  K,  p,  and  are  the  thermal  conductivity,  rlensity, 
and  heat  capacity  at  constant  pressure.  The  quantity  D - (K/ pC^)  is  also 
known  as  the  "diffusivity.  " For  almost  all  atmospheric  conditions  of  inter- 
est, heat  transfer  is  dominated  by  convection  rather  than  conduction.  This 
fact  is  equivalent  to  the  statement  that  negligible  energy  is  lost  by  a mole- 
cule in  the  time  required  for  the  molecule  to  move  across  the  optical  beam. 
Mathematically,  this  condition  can  be  expressed  as  V.  with  given 

by  eq.  (1).  Thus  for  a given  medium,  there  is  a minimum  value  for  the 
product  Vr^. 

To  ensure  that  the  blooming  distortions  5S  a function  of  propagation 
distance  are  the  same  in  the  laboratory  at  0.488  pm,  the  total  absorption, 
o-L,  and  the  Fresnel  m..mber,  must  be  held  constant.  The  absorption 
coefficient  is  o and  the  total  propagation  distance  in  the  blooming  medium  is 
L.  The  Frasnel  number  is  conveniently  defined  by 
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where  k = 2tt/\  and  \ is  the  wavelength. 

The  blooming  strength  is  measured'^  by  the  product  • P,p,  where 
P is  the  laser  power  at  the  entrance  to  the  blooming  medium  and  a.  is 
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defined  for  a gaseous  medium  by 
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The  quantities  y,  and  p are  the  molecular  polarizability  of  the  gas,  the 
specific  heat  (y  = C^/C^),  and  the  static  gas  pressure.  The  absorption 
coefficient  for  the  gas  is  a.  For  a liquid,  is 
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where  the  liquid  density  is  p,  the  nominal  index  of  refraction  is  n^,  and 
dn/dT  is  the  change  of  index  with  temperature. 

Wavelength  and  power  scaling  are  achieved  by  maintaining 
constant.  Table  3 presents  a scaling  example  comparing  10.6  pm  propaga- 
tion in  the  atmosphere  to  0.488  pm  propagation  in  a 50  cm  cell  filled  with 
Xe  and  NO^-  -As  can  be  seen  from  Table  3,  scaling  to  a 453  kW,  10.  6 pm 
blooming  experiment  can  be  achieved  within  the  laborator*^  using  reasonable 
experimental  parameters. 

As  mentioned  in  Section  II-C,  we  chose  xenon  seeded  with  NO^  for 
our  experiments  with  gaseous  blooming  media.  A partial  pressure  of  only 
3 to  4 Torr  of  NO^  is  required  to  achieve  50  to  70%  absorption  at  0.488  pm. 
The  cell  was  always  operated  with  a total  pressure  of  one  atmosphere.  We 
initially  felt  that  meaningful,  scalable  experiments  could  be  accomplished 
only  with  a gaseous  medium,  but  there  is  no  good  evidence  to  support  this 


Table  3.  Example  of  the  Scaling  Between  a 10.6  pm  Atmospheric 
Blooming  Experiment  and  a 0.488  pm  Laboratory  Gas 
Cell  Experiment 
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restriction  for  focused  beam  propagation.  As  data  in  the  next  section  show, 
similar  COAT  compensation  results  are  obtained  with  liquid  and  gaseous 
media  when  the  experimental  conditions  are  properly  scaled.  Sin  .e  strong 
blooming  is  much  easier  to  achieve  in  liquids,  there  is  little  motivation 
for  using  a gas  blooming  cell  for  future  experiments. 

2.  Experimental  Results 

Several  experimental  measurements  of  COAT  compensation 
for  blooming  have  been  made.  The  principal  differences  between  each  pair 
of  experiments  are  the  ratio  of  the  length  oi  the  blooming  medium  to  the 
focused  propagation  path  and  the  location  of  the  blooming  medium  within  the 
path. 

The  first  experiment  was  actually  part  of  the  Hughes  IR&D  program. 
The  experiment,  shown  schematically  in  Fig.  13,  used  a thin  liquid  cell 
filled  with  CCl^  doped  with  iodine.  There  was  no  transverse  wind  in  this 
experiment.  The  basic  ♦:wo-beam  arrangement  described  in  Section  II-E 
can  be  seen  in  Fig.  13.  When  the  liquid  cell  is  placed  as  close  as  possible 
to  lens  L^,  the  blooming  medium  is  in  the  collimated-beam,  near-field 
region  of  the  transmitting  array.  As  we  shall  see,  COAT  correction  has 
its  greatest  effect  on  near-field  blooming  distortions. 

With  the  cell  very  close  to  lens  L.  , the  results  shown  in  Fig.  14  are 
14  ^ 

obtained.  Both  the  boresight  irradiance  and  the  peak  irradiance  in  the 
beam  are  shown  as  a function  of  power  into  the  cell.  The  theoretical  curve 
shown  in  the  figure  is  obtained  from  a simplified  theory^' that  compares 
the  uncorrected  thermally-induced  beam  spreading  with  the  diffrac*-ion 
angle  of  a collimated  beam.  The  value  of  the  parameter  A is  obtained  by 
fitting  the  theory  to  the  experiment  at  one  power  level.  In  this  experiment, 
the  transmitted  beam,  is  not  focused  to  produce  the  far-field  irradiance  pat- 
tern, but  instead  expands  by  diffraction  to  produce  the  far-field  pattern. 

The  data  in  Fig.  14  indicate  that  the  COAT  system  increases  the 
optimum  transmitter  power’!'  by  a factor  of  2 to  3 and  increases  the  peak 
target  irradiance  by  up  to  a factor  of  4.  More  detailed  data  on  beam  pro- 
files  have  been  presented  elsewhere.  We  hypothesize  at  this  point  that  the 


’-The  transmitted  power  that  gives  maximum  target  irradiance. 
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-ig  14.  Thermal  blooming  compensation  with  a thin. 

blooming  cell  located  in  the  near-field  of  the  transmitter 
The  dasLd  curve  is  theoretical  (see  references  7 and  14). 
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COAT  system  is  this  effective  only  because  the  blooming  distortions  occur 
in  the  near-field  of  the  transmitter  and  that  as  the  blooming  occurs  farther 
from  the  transmitting  aperture,  the  COAT  phase  correction  will  be  less 

7 

effective.  (In  fact,  a simple  analysis  can  be  used  to  show  that  correction 
in  this  experiment  is  being  limited  mainly  by  the  finite  number  of  corrector 
elements  and  not  by  a fundamental  inability  of  COAT  to  remove  the 
distortions. ) 

To  test  this  hypothesis,  we  repeated  the  measurements  of  Fig.  14 

with  the  cell  at  different  distances  from  lens  L,  . The  results  are  shown  in 

2 

Fig.  15.  The  quantity  Zp^  = D.j,/4X.  is  the  classical  Rayleigh  distance, 
which  is  the  distance  at  which  the  transition  from  near-field  to  far-field 
occurs.  For  our  experiment,  Zp^  = 9 mm.  As  expected,  the  degree  of 
COAT  correction  is  reduced  significantly  as  the  blooming  cell  is  moved 
toward  the  far-field  of  the  transmitting  aperture  (d/Z„  >1).  In  fact,  it 

IV 

appears  that  significant  correction  occurs  only  for  distortions  located  well 
within  the  first  30%  of  the  propagation  distance  from  near-field  to  far-field 
since  no  appreciable  correction  is  observed  for  d/Zp^  > 0.28. 

The  next  experiment  employed  the  static  gas  cell  in  the  arrangement 
shown  in  Fig.  16.  In  this  case,  the  laser  beam  is  actually  focused  through 
the  cell.  The  data  obtained  with  this  arrangement  are  shown  in  Figs.  17 
and  18.  The  photographs  in  Fig.  17  are  black-and-white  photographs  from 

5 

a level-quantized  color  TV  display  with  1.  5 dB  per  color  level. 

There  are  five  principal  features  of  these  data.  First,  the  expected 
crescent-shaped  beam  distortion  is  observed  in  Fig.  17(c)  and  the  beam  is 
shifted  into  the  wind.  Second,  when  the  COAT  correction  is  applied,  the 
beam  maximum  is  maintained  on  the  boresight  axis  where  the  glint  is 
located.  Third,  there  is  no  improvement  in  peak  irradiance  when  the  COAT 
correction  is  applied,  except  for  a small  increase  in  the  power  through  an 
aperture  whose  width  equals  the  full-width  at  half-maximum  of  the  unbloomed 
beam.  Fourth,  the  system  convergence  time  was  unaffected  by  blooming 
level,  being  1.  5 to  2.0  msec  at  all  power  levels.  Finally,  the  convergence 
is  stable  at  all  power  levels. 
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Fig.  16.  Schematic  of  laboratory  optics.  Path  No.  1 through  the 

blooming  cell  and  path  No,  2 outside  the  cell  are  optically 
identical  when  the  cell  is  not  filled. 
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We  had  speculated  earlier  ' that  the  absence  of  correction  in  this 
experiment  may  have  been  caused  by  buoyancy  effects  or  by  too  low  wind 
velocities.  We  now  believe  that  we  see  no  correction  because  tlie  blooming 
medium  is  in  the  last  55%  of  the  propagation  path.  Strong  phase  distortions 
in  this  region  require  amplitude  changes,  not  phase  corrections,  at  the 
transmitting  aperture.  The  servo  response  speed  may  also  plfy  some  role 
(see  discussion  below),  but  this  is  not  clear  now.  However,  the  lack  of 
COAT  correction  in  data  in  Fig.  18  is  entirely  consistent  with  the  data  in 
Fig.  15. 

One  experimental  run  has  been  made  with  the  flowing  gas  cell.  The 
BAT  tracking  controls  were  included  in  these  tests.  The  experimental 
arrangement  was  similar  to  those  shown  in  Fig..  16,  except  that  lens  f^ 
(focusing  through  the  cell)  was  moved  and  its  focal  length  selected  so  that 
the  blooming  occurs  in  the  last  81%  of  the  propagation  path.  The  observed 
improvement  in  peak  irradiance  is  a factor  of  1.  35  or  35%.  This  improve- 
ment was  produced  by  the  phase  controls  alone;  the  tracking  controls  either 
alone  or  when  added  to  the  phase  controls,  produced  no  observable  change 
in  the  peak  target  irradiance. 

As  mentioned  previously,  we  were  not  able  to  determine  the  com- 
pensation of  thermal  blooming  using  the  dithered  focus  control  alone 
because  of  sphe’'ical  distortion  in  the  focus  actuator.  We  did  try  it  as  an 
addition  to  the  lli-channels  of  phase  control,  however,  but  we  observed  no 

change  in  the  peak  target  irradiance  when  the  focus  control  loop  was  closed. 

16  18 

This  result  is  consistent  with  other  work  that  has  demonstrated  little 

or  no  improvement  with  spherical  focus  control  alone,  but  significant  cor- 
rection with  coma  and  astigmatism  control. 

This  improvement  is  more  encouraging,  but  is  still  less  than  the 
factors  of  1.75  to  2.0  reported  by  others^^' The  comparison  is  unfair, 
however,  since  the  work  in  the  referenced  literature  used  slewing  beams 
(which  effectively  eliminates  blooming  near  the  focal  plane).  Our  experi- 
ments had  no  slewing,  and  the  distortion  was  in  the  last  81%  of  the  path, 
thus  presenting  a much  harder  task  for  the  COAT  system.  In  fact,  referring 
again  to  the  data  in  Fig.  15,  we  observed  a factor  of  1 . 7 improvement  in 
peak  target  irradiance  when  the  blooming  occurred  roughly  in  the  last  83% 
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of  the  propagation  path  (compare  the  uncorrected  curve  to  the  corrected 

curve,  d/ = 0.  17,  in  Fig.  15).  Actually,  in  the  experiment  leading  to 

Fig.  15,  the  blooming  occurred  only  from  d/ = 0.17tod/Zj^  = 0.28 

rather  than  over  the  last  83%  of  the  "focused"  path  (d/Z„  = 0.  17  to 

d/Z„  = 1.0).  Again,  this  is  a less  severe  case  as  far  as  COAT  correction 
K. 

is  concerned  than  the  experiment  giving  Fig.  19. 

The  final  blooming  experiment  used  a 20  cm  long  flowing  liquid  cell. 
Again,  an  optics  limitation  prevented  us  from  placing  the  entire  focused 
beam  path  within  the  blooming  medium.  The  experimental  arrangement  is 
similar  to  that  shown  in  Fig.  12.  Two  experiments  were  performed.  In 
the  first  experiment,  the  blooming  medium  was  in  the  last  72%  of  the 
propagation  path.  The  results  of  these  experiments  are  shown  in  Fig.  20. 
Not  unexpectedly,  the  data  are  very  similar  to  those  in  Fig.  19,  and  we 
observe  a factor  of  1.33  improvement  in  peak  target  irradiance.  This 
result  gives  us  confidence  that  there  are  no  significant  differences  between 
focused-beam  COAT  compensation  experiments  performed  with  liquid  and 
gaseous  media. 

In  the  secon'l  experiment,  the  liquid  cell  was  moved  to  the  first  72% 
of  the  path,  producing  the  data  in  Fig.  21.  We  now  see  a factor  of  1 . 5 
increase  in  the  optimum  peak  target  irradiance.  A factor  of  2 increase 
from  the  uncorrected  to  the  corrected  case  is  seen  at  the  highest  power 
level.  This  arrangement  produces  blooming  and  compensation  results 
along  the  propagation  path  that  are  similar  to  that  occurring  in  a slewed 
beam  when  the  entire  pre  agation  path  is  in  the  medium.^  ,19,20 
ments  using  stationary  and  slewed  beams  that  are  focused  through  the 
entire  blooming  medium  will  be  performed  as  part  of  follow-on  work  to 
this  contract  (contracts  F30602-76-C-0021  and  F30602-76-C-0022). 

3 . Conclusions 

We  now  feel  that  our  closed-loop,  real-time  phase  compensa 
tion  results  are  completely  consistent  with  other  results  that  were  obtained 
using  predictive,  open-loop  methods. A factor  of  about  1 . 5 to 
2.0  improvement  in  peak  target  irradiance  at  optimum  transmitter  power  is 
possible  with  multidither  COAT  dependong  on  the  slewing  rate.  Correction 
factors  of  2 or  more  may  be  possible  at  higher  transmitter  powers,  but  a 
lower  peak  target  irradiance  is  produced. 
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Fig.  19.  Thermal  blooming  compensation  data  obtained  with  flowing  gas 
cell.  Also  shown  is  the  focused  propagation  path  relative  to 
the  blooming  medium. 


454*  7R1 


8 cm  20  cm 

• — 4* H 


INPUT 


BLOOMING  t TARGET  PLANE 
MEDIUM 

(FLOWING  LIQUID) 


Fig.  20.  COAT  thermal  blooming  compensation  with  a flowing  liquid 

medium.  The  medium  is  in  the  last  72%  of  the  path.  Results 
are  similar  to  those  in  Fig.  19. 
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Fig,  21.  COAT  correction  with  a flowing  liquid  medium 
that  is  in  the  first  72%  of  the  propagation  path. 
Significantly  better  performance  is  observed 
than  in  the  experiment  described  in  Fig.  20. 


Care  must  be  exercised  when  comparing  compensation  results 

from  different  experiments,  since  our  data  indicate  that  only  the  blooming 

distortions  which  occur  in  the  first  20  to  30%  of  the  focused  propagation  path 

can  be  removed  by  correcting  the  transmitter  phase.  This  suggests  that 

any  method  that  Cc  n minimize  the  blooming  distortions  close  to  the  focal 

plane  will  be  effect  e not  only  in  increasing  the  target  irradiance,  but  also 

in  making  the  remaining  distortions  more  amenable  to  COAT  correction. 

Three  techniques  come  immediately  to  mind:  (1)  transmitter  intensity 

tailoring  or  contrcl^^  (2)  use  of  slewed  beams, (3)  use  of  pulsed 
25 

lasers  with  repetition  rates  low  enough  to  allow  only  2 to  3 pulses  during 
a wind-driven  medium  transit  time  across  the  beam.  When  comparing 
improvement  factors  in  different  experiments,  we  must  also  be  careful  to 
specify  whether  we  are  comparing  the  maximum  achievable  target  irradiance 
with  and  without  co”rection  or  the  improvement  at  a fixed  transmitter  power 
level.  The  use  of  slewing  beams  also  greatly  affects  the  comparisons. 

In  comparing  our  blooming  experiments  to  COAT  system  tests  done 

3 4 5 

with  no  distortions  or  with  atmospheric  turbulence  in  the  path,  w ■ have 

found  no  differences  in  system  convergence  time,  multiple  glint  discrimina- 
tion, or  tracking  ability.  We  checked  these  quantities  carefully  in  the  last 
two  experiments  described  above.  We  also  noted  that  the  same  blooming 
compensation  was  observed  with  a diffuse- surface  sphere  for  a target 
instead  of  a small  point-glint  target. 

Although  the  thermal  blooming  distortions  do  not  affect  the  JOAT 
servo  response  time,  we  are  concerned  about  what  the  response  time  should 
be  for  optimum  blooming  compensation.  Analysis  of  this  problem  is  greatly 
complicated  by  the  nonlinear  distortion  medium.  At  this  point,  let  us  just 
state  the  problem;  we  will  return  to  it  again  in  Section  IV. 

Since  the  COAT  system-plus-blooming  medium  is  a multiple-loop,, 
nonlinear  system,  it  is  not  unlikely  that  multiple  solutions  exist  for  the 
transmitter  phase  front  that  satisfies  the  servo  error  criterion:  near-zero 
dither  fundamental  frequency  in  the  signal  reflected  from  the  target.  The 
problem  is  to  choose  the  "path"  in  "correction  snace"  that  ensures  the  sys- 
tem stops  on  the  true  maximum  target  irradiance,  not  on  a "local"  maximum. 
It  is  important  to  point  out  at  this  juncture  that  we  are  only  hypothesizing 


the  existence  of  subsidiary  maxima  in  the  COAT  convergence  process  v/ith 
a nonlinear  distortion;  there  is  no  experimental  or  analytical  evidence  yet 
to  support  this  hypothesis. 

If  the  COAT  system  convergence  time  is  faster  than  the  response 
time  of  the  nonlinear  medium,  it  will  not  sample  the  true  slope  of  the  target- 
irradiance  versus  phase-error  curve  properly.  That  is,  the  instantaneous 
distortion  will  be  completely  removed.  If  instead  of  this  fast  response, 
the  COAT  system  was  slower  than  the  medium,  it  would  act  like  a "step, 
look,  decide,  step"  servo  system.  In  this  type  of  system,  a change  in  the 
transmitter  phase  front  io  made,  the  medium  responds  to  the  change,  and  a 
decision  is  made  whether  the  change  increased  or  decreased  the  target 
irradiance.  If  the  irradiance  increases,  the  direction  of  the  next  change  is 
the  same;  if  not,  the  direction  reverses.  This  is  effectively  the  type  of 
system  studied  by  others ^ The  amplitude  of  the  test  changes  (the 
dither  amplitude  in  a multidither  system)  will  also  be  important  in  this  type 
of  system  if  local  maxima  are  to  be  avoided. 

At  this  point,  we  speculate  that  it  may  be  important  to  allow  the 
medium  to  fully  respond  to  the  COAT-induced  phase  changes  and  possibly 
even  to  the  dither  frequencies  if  a multidither  COAT  system  is  to  converge 
to  the  global  optimum  transmitter  phase.  This  requirement  means  dither 
frequencies  well  below  the  present  8 to  32  kHz  band  used  in  the  DARPA/ 
RADC  system.  The  problem  with  such  a slow  system  is  that  it  is  totally^ 
inadequate  for  turbulence  compensation.  A possiole  solution  to  this  ' 
dilemma,  if  our  speculation  is  correct,  would  be  the  use  of  some  sort  of 
predictive  correction  so  that  the  COAT  system  starts  close  enough  to  the 
correct  phasefront  that  local,  sub-optimal  maxima  are  avoided.  A second 
solution  may  lie  in  the  existence  of  turbulence  in  the  optical  path.  The 
random  nature  of  turbulence,  and  a fast  COAT  system's  action  to  correct 
for  it,  may  be  enough  to  ensure  that  a broad  region  of  "correction  space" 
is  sampled  so  that  the  global  optimum  phase  correction  is  always  found. 

It  is  important  to  point  out  that,  all  of  the  above  discussion  aside, 
we  have  never  observed  any  servo  instability  in  our  experimental  bloom- 
ing compensation  measurements.  In  each  experiment,  the  servo  gain 
(response  time)  is  optimized  with  no  blooming  or  turbulence  distortions 


present.  The  distortions  are  then  added  and  the  laser  power  is  increased 
(for  blooming).  In  no  instance  did  the  COAT  servo  become  unstable. 

The  blooming  compensation  may  have  been  less  than  anticipated,  but  a stable 
phasefront  correction  was  always  achieved.  This  stab'  behavior  is  not 
unexpected  for  intensity-maximization  systems  such  xultiditner,  out- 

going wave  COAT,  but  may  not  be  characteristic  of  return-wave  systems 
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such  as  multidither,  return-wave  or  phase-conjugate  or  "TRIM-COAT"- 

systems. 

B.  Turbulence  Plus  Blooming  Compensation 

One  of  the  goals  of  this  contract,  and  an  important  consideration 
for  any  operational  COAT  system,  is  the  determination  of  whether  thermal 
blooming  will  diminish  a COAT  system's  ability  to  correct  for  turbulence. 

This  section  presents  the  results  of  our  studies  of  this  problem. 

1 . Artificial  Turbulence 

The  artificial  turbulence  phase  screen  was  produced  as  dis- 
cussed in  Section  II-D  and  located  as  shown  in  Fig.  12.  The  distortions 
proved  to  be  slightly  too  strong  for  this  COAT  system.  That  is,  more  and 
smaller  elements  are  required  to  fully  compensate  for  the  distortions. 

With  no  COAT  correction,  a time-averaged  (10  sec)  Strehl  ratio  of  0.22  was 
observed.  The  COAT  system  was  able  to  correct  this  beam  back  to  a Strehl 
ratio  of  0.  60.  The  phase  plate  was  rotated  to  produce  time-varying  turbu- 
lence. A plot  of  the  control  voltage  for  one  of  the  transmitting  array  ele- 
ments is  shovm  in  Fig.  22.  Significant  distortions  (as  measured  by  the 
corrections  for  them)  occur  at  frequencies  up  to  30  Hz.  This  is  in  reasonably 
good  agreement  with  the  50  Hz  spectral  width  observed  with  real  atmospheric 
turbulence  with  this  COAT  system. 

Using  a somewhat  less  severe  phase  screen,  we  have  compared  the 
use  of  phase  control  to  the  use  of  tracking  control  only.  This  phase  screen 
not  only  had  lower  phase  distortions,  but  the  quartz  plate  was  mounted  in 
such  a way  that  significantly  more  steering  error  was  introduced  than  with 

; TRIM-C0AT  is  the  acronym  applied  to  broad-band  image- compensation 
(I-COAT)  systems  that  are  used  to  correct  the  phase  of  a transmitted, 
coherent  beam. 
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the  phase  screen  shown  in  Fig.  12  and  used  in  the  turbulence  plus  blooming 
measurements  discussed  below. 

We  have  looked  at  four  cases:  (1)  no  COAT  control,  (2)  phase  con- 
trol only,  (3)  tracking  control  only,  (4)  phase  and  tracking  control  together. 
The  results  of  the  measurements  are  shown  in  Table  4.  As  can  be  seen 
from  the  data  in  the  table,  the  phase  and  steering  distortions  in  this  artifi- 
cial turbulence  screen  have  roughly  equal  contributions  to  the  uncompensated 
distortion. 

2 , Turbulence  Plus  Blooming 

The  turbulence  phase  screen  used  in  these  tests  was  the 
stronger  of  the  two  discussed  above.  The  thermal  blooming  scenario 
chosen  was  the  same  as  that  producing  the  data  in  Fig.  21:  a flowing  liquid 

medium  in  the  first  72%  of  the  propagation  path.  The  observed  compensa- 
tion performance  is  shown  in  Fig.  23.  The  peak  irradiance  for  turbulence 
plus  blooming  is  almost  as  good  as  that  achieved  in  the  absence  of  turbu- 
lence. The  difference  is  caused  by  the  inability  of  the  COAT  system  to  com- 
pletely remove  the  turbulence.  This  can  be  seen  by  noting  that  initial  slope 
of  the  ^ 0,  COAT-corrected  curve  in  Fig.  22  is  lower  than  the  slope  of 
the  free-space  line  and  by  noting  that  the  ratio  of  the  irradiances  for  the 
two  corrected  cases  at  any  power  level  is  very  close  to  0.60,  the  no-blooming, 
corrected  Strehl  ratio. 


Table  4.  COAT  Compensation  for  Artificial  Turbulence 


Case 

Peak  Target 
Irradiance, 
arb.  units 

Strehl  Ratio,  * 
S 

NO  COAT 

0.  39 

0.  65 

COAT  phase  control 

0.  52 

0.87 

COAT  tracking  control 

0.  52 

0.87 

Phase  and  tracking  control 

0.  57 

0.95 

=:<S  = 1 is  defined  for  no  turbulence,  COAT  phase  control. 
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Fig.  23.  COAT  correction  for  turbulence^n  the  presence 
of  blooming  distortions.  The  Cjyj  = 0 curves  are 
the  same  as  in  Fig.  21. 
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Conclusions 


The  data  presented  above  indicate  that  nearly  100%  of 
turbulence  distortions  can  be  removed  by  a COAT  System  (given  a sufficient 
number  of  COAT  channels)  even  though  the  absolute  Strehl  ratio  is  severely 
reduced  by  blooming.  The  nonlinear  blooming  medium  does  not  limit  the 
COAT  system's  ability  to  correct  for  turbulence.  1^..  conclusion  is  also 
supported  by  the  computer  simulation  data  discussed  in  Section  IV. 

C.  Complex  Target  Effects 

Contract  Amendment  No.  2 redirected  some  of  the  contract  work  to 
provide  time  and  funds  for  some  initial  studies  of  the  interaction  of  a multi- 
dither COAT  servo  system  with  backscatter  from  an  extended,  semidiffuse 

target.  Some  of  the  results  of  these  studies  were  presented  at  the  1975 

28 

Conference  on  Laser  Engineering  and  Applications.  The  interest  in 

pursuing  these  studies  was  initially  motivated  by  the  work  of  Ogrodnik  and 
29 

Gurski.  A discussion  of  the  results  of  these  studies  can  also  be  found 
in  Ref.  6. 

1 . Problem  Statement 

The  physical  problem  known  as  "speckle  noise"  can  be  stated 
briefly.  When  a coherent  laser  beam  illuminates  a target,  some  of  the 
energy  is  scattered  back  in  the  direction  of  the  transmitted  beam.  With  all 
except  single  glint  targets,  the  backscattered  radiation  produces  a random 
intensity  pattern  (a  "speckle  pattern")  at  a receiver,  which  is  located  close 
to  the  transmitter  for  cases  of  interest  to  us.  When  the  appropriate  condi- 
tions exist,  any  movement  of  this  speckle  pattern  relative  to  the  receiver 
will  produce  an  amplitude  modulation  of  the  received  signal.  Since  a multi- 
dither COAT  system  receiver  senses  amplitude  modulations  on  the  back- 
scattered  beam  that  are  produced  by  dithering  the  transmitted  beam  phase, 
large  amplitude  spurious  amplitude  modulations  in  the  receiver  may  swamp 
the  desired  modulations  and  thus  interfere  with  the  system  operation. 

Doppler  shifts  produced  by  rapidly  rotating  targets  may  also  produce 
false  signals,  but  these  spurious  signals  are  expected  to  be  less  important 
than  amplitude  modulation  effects  for  the  target  scenarios  of  interest.  This 


is  true  because  the  most  likely  glint  point  on  a rapidly  rotating  target  is  one 
nornaal  to  the  beana  incidence  direction;  this  point  will  have  nearly  zero 
transverse  doppler  shift.  Longitudinal  doppler  shifts  are  very  narrow  band 
and  thus  will  have  minimal  effect  on  a multiple- channel  COAT  system. 

Thus  for  COAT  applications,  't  is  important  to  determine  three 
things  about  target-signature  effects.  First,  what  kind  of  modulations  inter- 
fere with  the  proper  operation  of  a multidither  COAT  system?  The  interest- 
ing parameters  of  the  modulation  are  amplitude,  frequency  spectrum,  and 
power  spectral  density.  Second,  what  kind  of  receiver  modulations  are  pro- 
duced by  backscatter  from  real  targets  in  expected  operational  scenarios? 
That  is,  what  scenarios  will  produce  the  modulations  that  can  affect  the 
COAT  system?  Third,  can  a multidither  COAT  system  be  designed  that  is 
insensitive  to  speckle-induced  receiver  modulation  and  if  so,  what  are  the 
critical  design  parameters?  A follow-on  contract.  No.  F30602-76-C-0021 , 
to  the  current  contract  is  designed  to  produce  definitive  answers  to  all  these 
questions.  Some  preliminary  studies  of  the  magnitude  of  the  effects  were 
obtained  on  this  contract,  however. 

2.  Experimental  Results 

For  a first  experiment,  we  chose  the  simple  extended  moving 
target  shown  in  Fig.  24(a).  The  target  is  a cylindrical  scotchlite  strip  that 
could  be  rotated  at  different  rates.  The  target  is  larger  than  the  COAT- 
formed  array  pattern  in  both  dimensions  (parallel  with  and  perpendicular  to 
the  rotation  axis),  but  exceeds  the  element  pattern  diameter  only  in  one 
dimension.  The  curvature  of  the  scotchlite  provides  a localized  highlight'* 
on  which  the  COAT  system  can  lock.  The  laboratory  experimental  arrange- 
ment used  is  shown  in  Fig.  12.  No  turbulence  or  blooming  was  introduced 
in  the  paths  and  the  BAT  controls  were  not  used.  The  receiver /transmitter 
aperture  ratio  was  chosen  to  be  0.66. 


*It  is  erroneous  to  say  that  a multidither  COAT  system  requires  a "glint" 
for  convergence.  It  will  not  converge  on  a flat  plate  (at  least  the  systems 
considered  here  will  not),  but  any  highlight  a change  in  reflectivity  with 
beam  position  on  the  target)  will  suffice  to  provide  the  required  reference. 
We  have  established  this  fact  experimentally  with  several  target  shapes. 
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The  experiment  consisted  of  observing  the  peak  target  irradiance  and 
the  COAT  receiver  signal  frequency  spectrum  when  the  target  was  stationary 
and  moving.  Figures  24(b)  and  24(c)  shows  the  observed  COAT  photomulti- 
plier receiver  power  spectra  for  a typical  case.  The  spectrum  analyzer 
bandwidth  was  200  Hz.  Each  photograph  consists  of  five  2-min  scans  over 
a 0 to  50  kHz  range.  Figure  24(b)  is  the  observed  spectrum  when  the  target 
is  stationary  and  Fig.  24(c)  is  the  spectrum  when  the  target  rotates.  The 
target  rotation  speed  was  chosen  to  give  the  largest  observed  change  in  the 
spectrum  over  the  dither  band.  For  the  stationary  target,  the  noise  level 
is  low  and  the  dither  signals  are  clearly  evident  above  the  noise.  When  the 
target  is  moving,  the  noise  level  increases  significantly  and  most  of  the 
dither  signals  are  no  longer  evident.  The  peak  target  irradiance  was  the 
same,  however,  whether  the  target  was  stationary  or  moving.  For  this 
simple  target,  if  the  COAT  system  could  converge  on  the  target  when  it  was 
stationary,  the  system  could  converge  equally  well  when  the  target  was 
moving.  For  every  target  and  target  motion  rate  that  we  have  tried  to  date, 
this  same  observation  has  been  made:  no  noticeable  effects  when  the  target 
moves. 

The  observation  of  no  effect  on  the  COAT  system  performance  is 
contrary  to  a statement  made  in  the  third  management  report.  The  very 
earliest  experiments  did  seem  to  show  some  effect.  What  we  found,  how- 
ever, was  tliat  very  large  extraneous  dither  signals  were  being  introduced 
into  the  COAT  receiver  by  ground  loops,  not  by  target  speckle  effects. 

Since  the  ground  loops  have  been  eliminated,  we  have  not  been  able  to  find 
a target  which  "confused"  the  COAT  system  solely  because  of  its  motion. 

Since  we  have  seen  no  deleterious  COAT- system/ target-modulation 
effects,  a natural  question  is.  What  kind  of  spurious  signals  can  reduce  the 
system  performance?  We  thus  undertook  an  experiment  whose  goal  was  to 
generate  large  amplitude  spurious  amplitude -modulated  signals  capable  of 
ov^erloading  the  dither  servo  channels.  In  this  experiment,  the  target  was  a 
single  glint  and  a large  receiver  aperture  was  used.  Spurious  amplitude 
modulations  were  artificially  introduced  into  the  COAT  system  by  putting  an 
acousto-optic  modulator  (AOM)  in  the  optical  beam  path,  as  shown  in  Fig.  25. 
The  AOM  amplitude  modulates  the  transmitted  beam  at  a frequency  set  by 
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Fig.  25.  Experimental  setup  to  determine  effect  of  spurious  amplitude 
modulations  in  a COAT  receiver. 


the  voltage-controlled-os  illator  (VCO).  Modulation  depths  of  64%  were 
observed  over  the  dither  frequency  band  (8  to  32  kHz)  with  a triangular- 
wave  drive  to  the  VCO. 

Manually  tuning  the  VCO  through  the  dither  passband  sequentially 
overloaded  the  feedback  loop  at  each  dither  frequency  and  caused  a 10% 
reduction  in  peak  target  irradiance  (a  loss  of  proper  phasing  in  i of  the  18 
COAT  channels).  Spectrum  analyzer  measurements  indicated  that  the 
spurious  signal  level  of  100  mV  generated  in  this  manner  was  much  larger 
than  the  3 mV  dither  amplitudes.  This  initial  result  establishes  the  fact 
that  external  modulation  applied  to  the  COAT  beam  could  overload  the  con- 
trol circuitry  in  any  single  channel  and  cause  the  loss  of  that  channel  in 
forming  the  converged  beam. 

To  roughly  simulate  the  effect  of  worst-case  broadband  speckle 
modulations,  the  VCO  was  swept  over  a selectable  bandwidth  with  sweep 
periods  of  1 msec  and  0.  5 msec  (shorter  time  than  the  COAT  system  con- 
vergence time  of  1.5  to  2.0  msec) . During  this  series  of  tests,  the  COAT 
receiver  frequency  spectrum  was  recorded  and  the  power  level  in  the  con- 
verged beam  was  measured  with  the  sweep  voltage  turned  on  and  off.  The 
center  frequency,  f^,  was  selected  for  maximum  effect  on  the  COAT  system. 
A graph  of  the  number  of  active  elements  lost  versus  modulation  bandwidth 
is  shown  in  Fig.  26.  The  sweep  bandwidths  in  Fig.  26  range  from  1 kHz  to 
30  kHz  (f  = f^  ± BW/2),  but  the  ratio  of  the  peak  target  irradiance  with  the 

spurious  modulation  to  that  without  it  (P  ,/P  ) is  nearly  constant.  A 

2 mod  o 

value  of  P^^^/P^  = 0.  69  = (15/18)  corresponds  to  a loss  of  proper  phas- 
ing in  three  channels  of  the  average  and  (14/18)^  = 0.60  corresponds  to  a 
loss  of  4 channels.  Similar  results  are  obtained  when  the  sweep  period  is 
reduced  to  0. 5 nsec. 

The  power  spectrum  of  the  triangular-wave  swept  modulation  used 
here  consists  of  a sum  of  all  the  harmonics  of  the  sweep  frequency  spaced 
about  the  center  frequency: 


P = N a exp[i2TT(f  - nf  )]  , (4) 

m f j n ^ o sweep  •' 

ll  = -00 
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where  the  are  determined  in  part  by  the  sweep  bandwidth.  Although  this 

type  of  power  spectrum  will  probably  not  be  characteristic  of  any  target- 

induced  speckle  spectrum,  it  does  have  several  properties  which  a speckle- 

induced  receiver  spectrum  should  also  have.  First,  the  bandwidth  of  the 

extraneous  receiver  signals  is  limited,  but  will  almost  certainly  exceed 

29 

35  kHz  in  cases  of  most  interest.  Second,  the  product  of  peak  power 
spectral  amplitude  and  bandwidth  is  constant.  What  this  means  is  that  the 
receiver  can  see  a 100%  average  modulation-depth  signal  at  only  a single 
frequency;  as  the  modulation  bandwidth  increases,  the  average  receiver 
voltage  modulation  depth  (not  the  instantaneous  depth  or  amplitude)  goes 
down  as  the  reciprocal  of  the  square  root  of  the  bandwidth. 

Determination  of  the  actual  spectrum  of  a real  speckle-induced 
receiver  signal  is  beyond  the  scope  of  this  contract.  As  a worst  case, 
however,  consider  a spectrum  of  equal  amplitude  sinusoidal  frequencies, 
spaced  by  100  Hz  over  a 35  kHz  bandwidth,*'  which  is  close  to  a worst-case 

situation.  Since  there  are  350  frequencies,  the  average  receiver  voltage 

-1/2 

modulation  depth  for  any  one  must  be  (350)  w 0,  0^2,  or  4.  2%  to  main- 
tain a constant  powe r -bandwidth  product.  This  assumes  a 100%  speckle 
modulation  depth,  which  will  occur  only  for  purely  diffuse  targets.  If  we 
take  a 50%  modulation  depth  since  realistic  targets  will  probably  have  some 
specular  returns,  the  average  modulation  depth  drops  to  2,l%i.  This  value 
is  below  the  5,6%  peak  receiver  modulation  produced  by  a multidither  sys- 
tem operating  with  a ±20°  dither  modulation. 

A further  consideration  that  we  have  net  taken  into  account  is  that  the 
spectral  components  of  speckle-induced  receiver  noise  will  have  randomly 
varying  phases  with  respect  to  the  dither  modulation.  The  servo  synchro- 
nous detectors  will  thus  further  discriminate  against  this  type  of  noise.  We 
are  developing  an  analytical  demonstration  of  the  magnitude  of  this  discrimi- 
nation on  Contract  F3060Z-76-C-002  1 . Some  of  the  details  of  the  early 
stages  of  this  analysis  are  given  in  Rtg.  6 and  in  Section  IV-8  of  this  report. 


'■A  more  realistic  speckle-induced  spectrum  might  be  a uniform  triangular 
or  gaussian  distribution  with  a 50  kHz  bandwidth. 


3. 


Conclusions 


These  brief  studies  seem  to  indicate  that  target  - signature 
effects  are  not  going  to  be  as  severe  as  was  thought  earlier.  More  realistic 
targets  and  motions  need  to  be  studied,  however,  and  other  effects  such  as 
more  servo  channels  or  receiver  signal-to- shot  noise  need  to  be  investigated 
The  presence  of  turbulence  or  blooming  may  also  have  some  effect. 
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IV. 


ANALYTICAL  TECHNICAL  ACCOMPLISHMENTS 


Much  of  the  analytical  work  performed  on  this  contract  has  been 
accomplished  using  computer  simulation.  The  simulation  codes  model  both 
the  details  of  a COAT  servomechanism  (particularly  in  the  case  of  multi- 
dither COAT)  and  the  details  of  the  propagation  of  a coherent  beam  through 
turbulent  and  absorbing  medium.  The  codes  have  been  developed  to  their 
present  state  on  a variety  of  programs  including  the  Hughes  IR&D  program 
and  several  contracts:  N60921 -74-C-0249  (NSWC),  F30602-73-C-0248 
(DARPA/RADC),  F30602-75-C-0001  (DARPA/RADC).  The  results  of  the 
work  performed  on  this  contract  are  summarized  in  this  section. 


a 


A.  Computer  Simulation  of  COAT  Compensation 

1 • Summary  of  Capabilities  and  Limitations  of  the  Computer 

Simulation  Codes 

The  elements  of  the  adaptive  optics  computer  simulation  are 
indicated  in  the  flow  diagram  in  Fig.  27  The  simulation  contains  computer 
routines  that  model  the  operation  of  an  adaptive  phase  compensation  system, 
the  propagation  of  an  optical  beam  through  a turbulent  and  absorbing  medium, 
and  the  reflectivity  function  of  an  idealized  target.  The  adaptive  system 
simulator  consists  of  two  main  subblocks:  a servo  routine  that  can  model 
either  phase  conjugate  (return-wave)  control  or  multidither  (outgoing- wave) 
control  and  a phase  corrector  routine  that  can  model  segmented  or  continu- 
ously deformable  mirrors  as  well  as  tracking  and  focus  control.  The  latter 
capability,  added  under  Amendment  No.  1 to  this  contract,  is  discussed  in 
detail  in  Ref.  6.  The  simulation  models  a closed-loop  adaptive  system  in 
which  the  control  information  is  provided  by  the  detection  of  a return-wave 
reflected  from  the  target. 

The  effect  of  the  atmosphere  on  the  laser  beam  propagation  is 
modeled  by  a time-dependent  propagation  code  that  treats  the  effects  of  both 
turbulence  and  thermal  blooming.  The  time-dependent  code  is  used  on  both 
the  outgoing  and  return  paths.  The  difference  in  the  return  calculation  is 
that  the  propagation  is  modeled  by  a low-power,  linear  code,  which  includes 
the  effects  of  the  inhomogeneities  caused  by  turbulence  and  the  absorptive 
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Fig.  27,  Adaptive  optical  system  simulation  flow  diagram.  Both 
multidither  and  phase  conjugate  control  algorithms  are 
simulated  as  well  as  atmospheric  turbulence  and  thermal 


blooming . 


heating  introduced  by  the  high-power  transmitted  wave;  the  outgoing 
calculation  includes  nonlinear  heating  effects  (i.  e.  , the  outgoing  beam 
affects  the  medium). 

The  time-dependent  propagation  algorithm  used  in  the  simulation  is 
basically  a multipulse  code  in  which  the  laser  energy  is  propagated  in  a 
series  of  short  pulses.  The  pulse  length  is  assumed  short  enough  that  no 
heating  occurs  during  the  pulse.  Furthermore,  it  is  assumed  that  the 
heating  caused  by  a pulse  has  reached  a steady  state  by  the  time  the  next 
pulse  arrives.  The  term  steady  state  is  used  here  in  the  sense  that  all 
density  gradients  traveling  at  the  sound  velocity  are  assumed  to  have  left 
the  region  occupied  by  the  laser  beam  by  the  time  the  next  pulse  arrives. 
The  only  gradients  present  are  those  that  are  convected  with  the  local 
crosswind  velocity.  Each  pulse  "sees”  the  medium  created  by  the  steady 
state  heating  of  all  previous  pulses  but  is  unaffected  by  its  own  heating. 

This  model  allows  us  to  simulate  either  pulsed  or  cw  propagation  depending 
on  the  time  between  pulses  relative  to  the  time  it  takes  an  air  parcel  to 
traverse  the  beam.  Continuous  wave  propagation  conditions  are  simulated 
by  allowing  significant  overlap  between  the  re»t.t  ns  heated  by  successive 
pulses.  On  this  contract,  we  have  studies  only  cw  propagation. 

The  use  of  a time-dependent  propagation  code  differs  from  the 
approach  taken  by  Bradley  and  Herrmann  of  Lincoln  Laboratory,  who  have 
used  a cw  code  in  their  phase  compensation  studies.^  ^ Although  the  use  of 
a time-dependent  propagation  code  greatly  increases  the  complexity  and 
cost  of  the  calculation,  we  have  felt  compelled  to  use  this  approach  for 
two  reasons.  First,  the  time-dependent  approach  is  closer  to  the  physics 
of  the  striation  that  we  are  trying  to  simulate  since  we  are  primarily 
interested  in  adaptive  optical  systems  that  respond  in  a time  short  compared 
to  the  time  it  takes  the  atmosphere  to  change.  Such  adaptive  systems  per- 
form an  essentially  instantaneous  determination  of  the  required  phase  front. 
Moreover,  the  determination  is  continuously  updated  as  the  medium 
changes,  rather  than  to  be  put  in  and  held  fixed  until  a steady  state  is 
reached,  as  is  implicitly  assumed  when  a cw  propagation  code  is  used. 

The  distinction  between  the  time-dependent  and  cw  approaches  is  espe- 
cially important  for  the  simulation  of  adaptive  systems  that  employ 
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return-v/ave  algorithms?^  Whereas  it  can  be  argued  that  a cw  code  in 
conjunction  with  a maximization  routine  closely  models  the  behavior  of  an 
adaptive  optical  system  employing  an  outgoing-wave  multidither  algorithm, 
it  does  not  appear  possible  to  model  return-wave  systems  with  the  cw 
approach  since  such  systems  do  not  explicitly  maximize  the  irradiance  on 
the  target.  Instead,  they  continuously  set  the  phase  of  the  transmitted  wave 
equal  to  the  conjugate  of  the  phase  of  the  field  returned  from  the  target. 
There  does  not  appear  to  be  any  reason  to  expect  that  the  final  state  reached 
by  such  a system  can  be  obtained  from  a cw  calculation  or  a sequence  of 
such  calculations. 

The  other  reason  for  the  use  of  a time-dependent  propagation  code 
in  our  studies  is  the  desire  to  follow  the  time  development  of  the  correction 
process  especially  in  the  case  of  adaptive  systems  employing  outgoing- wave 
multidither  control  algorithms.  An  outgoing- wave  multidither  control  sys- 
tem attempts  to  maximize  the  irradiance  at  the  target  at  each  instant.  Vv' e 
wish  to  know  if  this  instantaneous  maximum  control  policy  leads  to  a global 
maximum  in  the  thermal  blooming  correction  problem  or,  instead,  yields 
a secondary  maximum.  This  can  be  determined  by  comparing  the  maximum 
obtained  with  the  time-dependent  and  the  cw  codes  because  the  cw  approach 
yields  the  global  maximum. 

The  time  sequence  of  events  in  the  simulation  is  shown  in  Fig.  28. 

A high-power  pulse  of  laser  energy  is  propagated  to  the  target  and  the 
resulting  change  in  the  density  of  the  air  in  the  region  occupied  by  the  beam 
is  calculated,  taking  into  account  the  convective  effect  of  the  local  crosswind 
velocity.  It  is  assumed  that  the  adaptive  phase  compensation  is  applied  just 
prior  to  the  firing  of  the  next  high-power  laser  pulse.  Moreover,  we  are 
presently  assuming  that  the  adaptive  optical  system  has  infinite  bandwidth 
so  that  it  is  able  to  determine  and  apply  the  phase  correction  instantaneously 
without  any  delay  oi  Lime-averaging  effects.  The  information  required  to 
determine  the  phase  correction  is  obtained  by  transmitting  a low-power 
reference  pulse  to  the  target  and  reflecting  the  same  from  the  target.  This 
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Fig.  28.  Time  sequence  of  events  in  multipulse  adaptive 
optics  simulation: 

(1)  No  heating  during  a pulse  (t^  <<  t^^oustic^ 

(2)  Steady  state  heating  for  t-t^  > 

(3)  Adaptive  optics  correction  applied  in  the 
interval 


reference  pulse  propagates  through  the  medium  that  will  be  seen  by  the 
next  high-power  pulse  and  thus  provides  the  necessary  phase  correction 
information.  This  process  is  then  repeated  for  each  subsequent  high-power 
pulse. 

With  the  exception  of  some  runs  in  which  the  effect  of  the  atmosphere 
was  represented  by  a single  phase  screen  in  front  of  the  ti  :>nsmitter , all  of 
the  data  that  have  obtained  to  this  date  apply  to  the  simulation  of  adaptive 
systems  utilizing  a return-wave  control  algorithm,  i.  e.  , systems  that  cor- 
rect by  inserting  the  conjugate  of  the  phase  of  the  wave  returned  from  the 
target.  The  reason  for  this  concentration  on  the  phase  conjugate  simulation 
is  simply  that  such  simulations  are  significantly  easier  to  ir  rplement  and 
less  costly  to  run  than  outgoing- wave  simulations.  The  simulation  of  a 
return-wave  system  requires  simply  that  we  transmit  a sit  gle  low-power 
reference  wave  to  the  target  and  a single  low-power  reflected  wave  back  to 
the  transmitter  aperture.  The  phase  correction  is  then  equated  to  the  con- 
jugate of  the  phase  of  the  reflected  field.  In  contrast,  an  accurate  multi- 
dither simulation  requires  the  sending  and  receiving  of  many  low-power 
reference  pulses  for  each  high-power  pulse.  We  must  sample  the  modula- 
tions of  the  integrated  irradiance  of  the  target  returns  at  a rate  set  by  the 
highest  dither  frequency.  For  example,  it  was  found  in  our  earlier  multi- 
dither servo  system  simulations  that  a sampling  time  increment  of  5 psec 
is  required  for  an  18-channel  system  having  an  upper  dither  frequency  of 
3?,  kHz  and  an  overall  convergence  time  of  1 msec.  Hence,  approximately 
20t  samples/ convergence  time  were  required  in  those  simulations.  It  was 
practical  to  i -ike  200  samples / convergence  lime  in  the  earlier  work  because 
of  the  simplicity  of  the  propagation  model  that  was  used.  In  the  present 
case,  however,  the  propagation  model  is  quite  complicated  and  relatively 
costly  to  run  and  thus  it  is  clearly  not  practical  to  do  an  equivalent  simula- 
tion. Additional  details  on  the  computer  simulation  codes  — servo,  propaga- 
tion, mirror,  and  target  — are  given  in  Ref.  30. 

2.  Turbulence  and  Blooming  Compensation 

30 

We  have  demonstrated  elsewhere  that  the  effect  of  atmo- 
spheric turbulence  on  the  ensemble-averaged  irradiance  distribution  can  be 
accounted  for  by  introducing  a random  phase  screen  in  front  of  the  transmitting 
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apertur  i.  The  appropriate  phase  distribution  is  that  which  would  be 

'\C 

produced  by  a point  source  located  on  the  target.'  This  implies  that  an 

ideal  adaptive  optics  system  should  be  able  to  completely  compensate  for 

the  effects  of  the  turbulence  by  introducing  the  conjugate  phase  distribution 

as  a correction.^  The  degree  to  which  this  is  achieved  depends  on  whether 

the  target  has  a point-glint  scatterer  that  provides  the  requisite  phase 

information  and  on  the  degree  to  which  the  adaptive  optics  system  can 

reproduce  the  desired  phase  distribution. 

30 

Previously  we  have  demonstrated  that  phase  conjugate  COAT  con- 
trol with  a deformable  mirror  can  correct  a beam  focused  through  a turbu- 
lent atmosphere  back  to  within  70  to  95%  of  the  diffraction-limited  (free- 

30 

space)  target  irradiance.  We  have  also  demonstrated  analytically  that  a 

multidither  COAT  system  will  arrive  at  the  same  phase  correction  as 

a phase  conjugate  system  when  correcting  for  linear  (not  power-dependent) 

distortions.  The  amount  of  correction  depends  on  the  range  (Z),  turbulence 

strength  (C^),  wavelength  (\),  transmitter  diameter  (D^),  and  number  of 

deformable  corrector  mirror  actuators  elements  (N  ).  The  dependence  of 

^ 6 30 

the  residual  phase  ervor  after  correction  is  described  by  ' 

= 0.05l(^)^  C^  Z (5) 

where  <4>g>  is  the  residual  mean  square  phase  error  acrosf  the  aperture 
after  correction.  This  result  includes  an  approximation  for  the  shape  of  a 
deformable  mirror  surface.  The  effective  Strehl  ratio  is  given  by 

S a exp  • (6) 

2 

The  value  of  S is  plotted  in  Fig.  29  versus  Cj^  for  several  cases  of  interest. 


-‘Tn  this  discussion,  we  assuixie  that  amplitude  scintillation  effects  are  small. 
If  they  are  not  small,  the  effect  of  atmospheric  turbulence  cannot  be  repre- 
sented simply  by  a phase  screen  but  rather  we  must  also  introduce  a ran- 
dom apodization. 

+ 

'Subject  to  the  assumption  that  amplitude  effects  are  negligible. 
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Fig.  29.  Strehl  ratio  associated  with  atmospheric  turbulence  ver 
BUS  number  of  mirror  actuators. 
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To  provide  a check  on  the  above  analytical  work,  we  have  perforrned_ 
two  types  of  mirror  simulations.  In  each,  the  effect  of  a turbulent  phase 
screen  placed  in  front  of  an  aperture  was  compensated  by  a mirror  with  a 
finite  number  of  actuators.  The  average  Strehl  ratio  obtained  from  a 
sequence  of  five  independent  phase  screens  was  determined  as  a function  of 
the  number  of  actuators.  In  one  of  the  simulations  the  deformable  mirror 
software  was  used  in  conjunction  with  a sinusoidal  multidither  COAT  algo- 
rithm of  the  type  used  previously  at  HRL  in  our  COAT  servo  system  studies. 
In  the  other  simulation,  the  mirror  was  modeled  by  a segmented  mirror  with 
piston  and  tilt  control  on  each  segment.  The  piston  and  tilt  settings  for  each 
segment  were  determined  by  a least  squares  fit  to  the  random  phase  surface 
over  the  segment. 

The  agreement  between  the  piston  and  tilt  results  and  the  theoretical 
predictions  obtained  from  eqs.  (5)  and  (6)  is  very  good  (see  Ref.  30),  which 
gives  us  confidence  in  the  theoretical  results  given  in  these  equations.  We 
believe  that  these  results  are  representative  of  those  that  would  be  obtained 
with  a deformable  mirror  in  the  absence  of  the  2Ntt  problem,  is  avoided. 

The  agreement  between  tjie  theory  and  the  deformable  mirror- 

multidither  COAT  simulation  results  is  reasonably  good  for  moderate  values 

2 2 

of  Cj^Z  but  is  poor  for  large  values  of  C^^Z.  We  attribute  the  poor  results 
obtained  at  large  values  of  C^Z  to  a 2NTT-type  of  behavior.  A similar  prob- 
lem is  observed  with  the  focus  control  computer  simulation.  The  deform- 
able mirror  simulation  that  was  used  in  these  runs  has  a 2Nit  correction 
loop  that  introduces  a 2tt  correction  whenever  the  phase  difference  between 
actuators  exceeds  four  radians.  The  intent  is  to  suppress  2tt  errors  intro- 
duced by  the  servo  system.  However,  if  the  phase  distortion  that  is  to  be 
removed  by  the  mirror  changes  by  more  than  four  radians  between  actua- 
tors, the  2it  "correction"  introduced  by  the  2Ntt  correction  loop  is,  in  fact, 
a 2tt  error  and  the  m.irror  performance  is  corresponding  degraded.  This 
problem  could  be  avoided  by  removing  the  2Ntt  correction  loop  but  then  we 
would  be  faced  with  2Nit  servo  errors.  One  way  to  avoid  this  problem  is 
to  use  more  actuators  so  that  the  phase  change  between  actuators  never 
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exceeds  four  radians.  Another  way  would  be  to  design  the  mirror  so  that 
the  likelihood  of  2Nit  errors  is  reduced,  in  which  case  the  2Nti  correction 
loop  could  be  eliminated.  More  work  clearly  needs  to  be  done  on  this 
problem. 

V/e  have  made  several  comp\:.ter  runs  to  determine  turbulence  com- 
pensation using  only  the  tracking  and  focus  controls.  Figure  30  shows  the 

observed  performance  with  tilt  control  alone.  The  relationship  between  the 

2 

quantity  P in  Fig.  30  and  Cj^Z  is  shown  in  Fig.  31  for  several  cases  of 
interest.  Figure  32  shows  a summary  of  the  compensation  results  achieved 
with  BAT  controls  compared  with  those  produced  by  a 37-element  multi- 
dither COAT  control  system.  The  deformable  mirror  used  in  all  the  simu- 
lation runs  has  37  actuators.  The  assumed  propagation  scenario  is  also 

2 - 1 1 1/3 

shown  in  the  figure.  For  reference,  S 21  x 10  m i j considered 

to  be  very  strong  turbulence.  Also  shown  in  the  figure  is  the  theoredcal 
limit  given  by  eqs,  (5)  and  (6)  for  a 37-element,  deformable  mirror,  phase- 
conjugate  COAT  system  that  has  perfect  phase- sensing,  but  which  is  limited 
by  the  ability  of  the  mirror  to  provide  th.;  desired  phase  front. 

There  are  several  noteworthy  features  of  the  data  in  Fig.  32.  First, 
the  agreement  between  analysis  and  computer  simulation  for  the  37-element 

deformable  mirror  COAT  system  is  very  good.  The  slight  discrepancy  at 
2 

large  values  of  0^^  Z may  be  caused  by  improper  operation  of  the  simulation 
routine  designed  to  remove  2mr  ambiguities,  as  discussed  above.  Second, 
significant  correction  is  obtained  vdth  the  tracking  or  focus  controls  when 
used  separately,  but  the  two-axis  focus  control  provides  the  most  correc- 
tion. In  fact,  when  the  tracking  controls  are  combined  with  the  two-axis 
focus  control,  little  increase  in  Strehl  ratio  is  produced  over  that  observed 
with  the  focus  control  alone.  We  attribute  the  lack  of  a significant  tilt  com- 
ponent in  the  simulation  (and  thus  little  correction  with  the  tracking  control) 
to  our  assumption  of  an  infinite  outer  scale  for  the  turbulence.  If  we  had 
used  a finite  outer  scale,  we  expect  that  the  tracking  correction  would  have 
been  more  significant. 

The  significantly  worse  performance  of  a spherical"''  focus  control 
2 - 1 1 1/3 

above  Cp^^Z  = 40  x 10  m is  also  interesting.  Finally,  not  unexpectedly, 

'i=Two-axis  focus  control  with  both  axes  driven  by  the  same  dither  and  cor- 
rection signals. 
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a 37-channel  multidither  COAT  system  produces  significantly  more 
correction  than  the  BAT  controls  alone;  the  BAT  controls  produce  roughly 
factor  of  2 improvement  in  Strehl  ratio  and  the  37-channel  COAT  system 
can  improve  the  Strehl  ratio  another  factor  of  2. 


8 


a.  Turbulence  and  Blooming  — As  we  reported  earlier, 

the  effects  of  turbulence  and  thermal  blooming  can  be  considered  separately, 
and  then  combined  in  a mean-square  sum  fashion  obtain  the  beam  Strehl 
ratio  when  both  effects  are  considered  simultaneously.  With  this  result,  it 
is  not  too  surprising  that  COAT  compensation  for  turbulence  is  unaffected 
by  the  presence  of  thermal  blooming.  The  experimental  data  presented  in 
Section  III  demonstrate  this  conclusion.  Our  computer  simulations  also  con- 
firm that  good  turbulence  compensation  is  possible  even  if  little  or  no  bloom- 
ing compensation  occurs.  Two  examples  of  the  simulation  results  are  shown 
in  Figs.  33  and  34.  The  transmitted  beam  in  this  case  is  a gaussian, 
truncated  at  the  10%  intensity  radius. 

There  are  two  principal  observations  to  be  made  on  the  data  in 
Figs.  33  and  34.  First,  there  is  very  little  thermal  blooixting  compensation. 
This  result  is  not  consistent  with  our  experimental  work  nor  with  other 
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experimental  ' or  computer  simulation  work.  We  feel  that  the 


absence  of  correction  in  our  simulations  is  caused  by  (1)  the  use  of  a return- 


vave  correction  algorithm  and  (2)  by  an  isoplanatic  problem  associated  with 
30 

unite-width  glints.  Our  experimental  data  were  taken  with  a fast  outgoing- 
ve  multidither  system.  The  data  obtained  by  others^^’  I?,  19 


v/a 


produced  by  outgoing- wave  systems,  but  in  each  case,  the  system  was 
effectively  a slow  dither  system. 

We  are  beginning  to  accumulate  evidence  that  the  inability  of  a phase 
conjugate  system  to  correct  for  blooming  — and  its  tendency  to  decrease  the 
target  irradiance  with  strong  blooming  - may  be  common  to  all  "return- 
wave"  COAT  systems.'''  A return-wave  COAT  system  is  one  that  operates 
only  on  the  wavefront  returned  from  the  target  to  derive  the  servo 


■'We  have  observed,  as  have  others,  that  for  transmitter  powers  exceed- 
ing the  optimum  power  for  blooming,  a phase  conjugate  COAT  correction 
is  worse  than  doing  nothing  at  all. 
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33.  Phase  conjugate  COAT  correction  for  thermal  blooming  at 
10.6  ^Jlm  for  a truncated  gaussian  beam  with  beam  slewing. 
Cases  with  and  without  turbulence  are  shown. 


3719-«3n2 


Fig.  34.  Phase  conjugate  COAT  correction  for  thermal 

blooming  and  turbulence  at  3.8  [xm  for  a truncated 
gaussian  beam  with  beam  slewing. 
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correction  signals.  Phase-conjugate  systems,  return-wave  mull  idither 
systems,'  and  hybrid  "TRIM-COAT"  systems*"  are  all  reUirn-wave  systems. 
In  fact,  of  the  presently  known  adaptive  systems,  only  the  multidither, 
outgoing- wave  COAT  system  appears  to  be  a true  maximization  system. 

Our  simulation  work  has  also  indicated  that  referencing  off  of 
extended  target  regions  in  the  presence  of  blooming  will  give  incorrect 
correction  results.  The  problem  is  the  thermal  blooming  equivalent  of  the 
finite  isoplanatic  patch  size  encountered  in  trying  to  image  extended  targets 
through  turbulence.  If  the  target  reference  region  exceeds  an  atmospheric 
coherence  diameter,t  a COAT  system  will  not  be  able  to  obtain  the  correct 
phase  errors  for  the  focused  propagation  path.  We  have  observed  during 
our  computer  simulation  work  with  phase  conjugate  correction  for  blooming 
that  this  problem  occurs  even  with  glints  less  than  one-half  the  diameter  of 
the  diffraction-limited  transmitted  beam  on  the  target.  We  plan  to  study 
this  problem  further  on  other  related  programs. 

The  second  observation  on  the  data  in  Figs.  33  and  34  is  that  the 
ability  to  compensate  for  the  effects  of  turbulence  is  apparently  not  degraded 
by  the  presence  of  thermal  blooming  effects.  The  net  improvement  in  beam 
quality  with  phase  conjugate  compensation  is  most  impressive  for  the  3.8  pm 
case  in  Fig.  34  and  for  cases  even  at  10.6  pm,  with  large  transmitter  aper- 
tures. This  is  to  be  expected,  of  course,  since  turbulence  effects  are  most 
pronounced  in  the  cases  of  large  beams  and  shorter  wavelengths. 

B.  Complex  Targets  and  Speckle  Modulations 

As  part  of  the  contract  redirection  defined  by  Amendment  No.  2, 
we  have  been  pursuing  a study  of  the  behavior  of  a multidither  COAT  system 
when  it  is  subjected  to  multiplicative  noise  in  the  receiver.  The  type  of 


*The  reflected  coherent  light  sensed  by  the  receiver  is  dithered,  rather 
than  the  transmitted  beam. 

Transmitting -Imaging;  broad  band,  incoherent  target  returns  are  sensed 
by  an  image  - compensating  Imaging-COAT  system  and  this  information  is 
used  to  correct  the  phase  of  the  transmitted  beam. 

tThe  distance  between  two  points  in  the  observation  plane  for  which  rays 
emanating  from  a point  source  see  the  same  optical  path  length. 
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noise  considered  is  very  similar  to  the  modulations  produced  by  complex- 
target  speckle  effects  (see  Section  III  and  Ref.  Z9).  Our  analytical  studies 
have  utilized  both  closed  form  analysis  and  computer  simulation  of  an 
18-channel  multidither  servo. 

1 . Analytical  Models  for  Speckle  and  COAT  Servo 

In  the  previous  quarterly  report,  ^ we  developed  a statistical 
analysis  intended  to  predict  and  explain  the  behavior  of  a COAT  system  sub- 
jected to  multiplicative  speckle  modulations.  This  analysis  will  now  be 
used  to  calculate  expectation  values  of  various  COAT  system  parameters. 
We  then  compare  these  calculations  to  data  gathered  from  computer 
simulations . 

The  first  parameters  of  interest  are  P^,  the  average  power  of  the 
dither  modulation  signal,  and  Pg,  the  average  power  of  the  speckle  signal 
within  ±400  Hz  of  any  dither  frequency.  The  expression  derived  in  Ref.  6 
for  PjQ  for  N servo  channels  is 


n 


f 


where  K is  a receiver  proportionality  constant,  is  the  dither  amplitude, 
and 


N 

'■t.  = ^ -Pm>  • 

m=l 
n^m 

The  J (v)  and  J,(4>)  functions  in  eq.  (7)  are  the  zero-  and  first-order  Bessel 

® ^ th 

functions  and  is  the  phase  error  signal  in  the  i control  channel.  The 

expectation  value  of  Pj^  is 
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n-^1 


(8) 


Since  is  a function  of  the  random  variables  Pj , p^*  • • • # it  is  also 

a random  variable  whose  expectation  value  is  expressed  as 


<b 


dp^.P(Pi)P(P2)--  P(Pj,)  b^(Pi.  P^---  Pi,) 


(9) 


where  P(P.)  is  the  probability  density  of  p^. 

We  assume  that  the  p.'s  are  uniformly  distributed  between  the  phase 

angles  a and  a + a as  indicated  in  Fig.  35.  Since  a is  arbitrary,  it  is 
“ o o o 

convenient  to  define  it  as  a = 'v/2;  this  will  make  the  integrals  in  eq.  (9) 
easier  to  calculate.  Physically,  a is  the  instantaneous  mean  phase  error, 
which  the  COAT  system  is  trying  to  drive  to  zero.  With  this  definition,  as 
a approaches  zero,  the  COAT  array  converges.  Conversely,  as  a approaches 
2tt,  the  array  approaches  a convergence  level  of  zero. 

For  any  given  p^,  the  probability  density  is  P(Pj)  = over  a region 

-a/2  < < a/2,  and  zero  elsewhere.  Using  the  above  information  plus  the 

definition  in  eq.  (7),  we  can  write  eq.  (9) 


<b^  = a-’ 
n 


N N ralz  roll  roll 

•‘^K  -Pm>  -Pk' 

m = l k=l  J-o!Z  J-alZ  J-qIZ 


m = l k=l  J^or/ 2 
m^k,  n 
k^n 


dp  sin 


n 


ip„-p 


m 


N rafZ  ro/Z 


(10) 


Fig.  35.  Phaser  diagram  of  a partially  con- 
verged, N-element  COAT  array 
defining  the  quantities  o and  used 
in  the  text. 


Carrying  out  the  integrals  and  summations  in  eq.  (10)  gives 


<b^>  = i (N^  - 3N  + 2) 
n z 


L‘-  » J 


1 


[i . 


+ i<N  - 1) ^ 


(11) 


Substitution  of  eq.  (11)  into  eq.  (8)  gives  the  expectation  value  of  the  average 
dither  p 'wer: 


The  expectation  value  of  P will  now  be  calculated.  As  derived  pre- 

1 6 ® 
viously , 


P 

s 


N M 


k=l  j=l 


(13) 


where  I is  the  mean  convergence  level  and  a.k  is  a random  amplitude 
m J 

variable  that  characterizes  the  amplitude  of  the  spurious  modulation  at  fre- 


quency OJ 


jk' 


We  have  assumed  that  there  are  M discrete  modulation  fre- 


quencies around  each  of  the  N dither  frequencies  so  that  the  spurious  modula- 
tion has  the  form 


(14) 


M 


N 

I ^ZlLi'jk'v^v 

k-l  j=l 


where  4>.k  is  a randcmly  chosen  phase  between  zero  and  Ztt. 


,2  . 


Since  Ij^  is  a function  of  the  independent  random  variables 
• • • Pjvj  ^ke  a^j^  are  themselves  independent  random  variables,  the 


expectation  value  of  can  be  expressed  as 


<P  > 
s 


i <l2 


N M 
k=l  j=l 


(15) 


In  our  computer  simulations,  the  a^j^  were  chosen  randomly  between  the 


values  of  0 and  a 


Tne  probability  distribution  is  this  uniform  over  that 
interval  and  zero  elsewhere.  The  expectation  value  of  their  square  is 


2 

<a.,  > 
Jk 


max 


2, 

a da 


max 


(16) 


max 


The  meein  convergence  level,  Ij^,  is  calculated”  to  be 


I,,  = 


t2...  N N 

J (4^) 


"M 


N _2 


N 


2^  2^  cos(p^-p 
n=  I m=  1 
n^m 


m 


(17) 
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The  expectation  value  of  Ij^  is 


/ipi  • /p 


dp^  P(Pl)  P(P2>  • • • P(Pn>  Im <Pi-  Pz  • • • Pn> 


Substituting  the  definition  of  given  in  eq.  (17)  and  the  probability  densities 
(P(p^)  = d ) into  eq.  (18),  we  can  solve  the  expectation  value  of  Ij^.  These 
steps  will  be  omitted,  however,  for  brevity.  The  result  of  the  calculation  is 


. 2 


r V N / (a/ 2) 


/.  6 11  6 \ sin  (q/Z) 

A N (»/2)^ 


+ 2 4- 


(?-?)[' -“S"]}  ■ 


(19) 


Figure  36  plots  the  root-mean- square  value  of  found  from  eq.  (19)  as  a 
fraction  of  a for  an  18-chcinnel  COAT  system.  When  eqs,  (16)  and  (19)  are 
substituted  into  eq.  (15),  the  expectation  value  of  the  average  power  in  the 
speckle  modulation  signal  is  found  to  be 


<?s>  = -g-  K^(NM) 


max 


1 

'n^ 
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+ 2J  (V(j)  f -f.  T - T 


(”  ■ n") 


sin 


(0/2) 


(o/Z)‘ 
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[(' 


k.  + li_  6 \ sin"^  (p/2) 

^ N^'nV  (W2)4 


(20) 


Expectation  values  of  many  other  parameters  can  be  calculated  in  a 
similar  manner.  Each  calculation  will  give  expressions  that  are  functions  of 


a = CONVERGENCE  ANGLE 


Statistical  expectation  value  of  the 
mean  convergence  level,  for  an 

I 8-channel  multidither  COAT  system 
(dither  amplitude  = 20  ),  calculated 
from  eq  . (13). 


N,  lij,  and  o and  possibly  other  variables  such  as  a , K.,  etc.  , which  are 
included  in  the  definitions  of  a given  parameter.  For  any  particular  COAT 
system,  most  of  these  variables  will  be  specified,  The  variable  a,  how- 
ever, is  a statistical  parameter  which  characterizes  the  convergence 
level  based  on  an  intuitive  model.  It  is  therefore  questionable  whether  one 
could  ever  measure  such  a parameter  experimentally.  Nevertheless,  know- 
ing several  expectation  values  as  functions  of  allov/s  us  to  derive  expres- 
sions for  one  expectation  value  as  a function  of  another:  <Pq>  as  a function 
of  <Ij^>»  example.  This  kind  of  analytical  prediction  can  be  checked 
against  the  real  world  and  is  therefore  more  useful. 

Figure  37  shows  plots  of  <Pg>  and  versus  <Ij^>  calculated  in 

this  manner.  Several  curves  of<P^>  corresponding  to  different  values  of 

b 

a are  shown.  One  can  appreciate  how  speckle  modulations  effect  a 
max 

COAT  system  by  consideiing  one  value  of  a , e.  g, , = 0.  006. 

Below  a convergence  level  of  approxlrriately  0.  5 the  power  in  the  dither 

signal  is  greater  than  the  power  in  the  speckle  signal  while  the  converse  is 
true  above  this  level.  Based  on  the  simulation  data,  the  COAT  system 
apparently  has  no  trouble  converging  to  = 0,  S,  where  power  ratio, 
p = <Pq>/<P  >,  is  equal  to  unity.  However,  when  a convergence  level  of 
approximately  1.  77  is  reached,  the  convergence  process  stops  and  a random 
oscillation  between  the  convergence  level  of  0.  5 and  0.  7 7 begins  (the  com- 
puter simulation  runs  discussed  below  exhibit  this  behavior). 

From  the  data  in  Fig.  37  it  is  evident  that  for  any  given  a the 

m3>X! 

ratio  Rj,  of  the  dither  power  to  the  speckle  power  varies  over  a large 

range  of  values  as  a function  of  convergence  level.  For  the  example  of 

a = 0.  006.  the  power  ratio  is  very  small  at  a convergence  level  of 
max 

O.  77  (the  actual  value  is  about  0.  3).  With  the  speckle  power  so  much 
larger  than  the  dither  power  it  is  reasonable  to  expect  convergence 
problems. 

2 . Computer  Simulation  Code 

Along  with  this  cinalysis,  we  have  been  making  computer 
runs  using  the  detailed  COAT  servo  simulation  shown  in  Fig.  3 8.  The 
spurious  modulation  function  (as  defined  by  eq.  (14))  in  our  work  to  date) 


is  introduced  as  multiplying  the  normal  receiver  signal.  For  the  simulation 

studies  performed  to  date,  we  have  chosen  M = 21  and  spaced  these  noise 

frequencies  ±83  Hz  around  each  dither  frequency.  This  is  a rather  arbitrary 

choice,  but  it  does  ensure  that  all  spurious  modulation  is  well  within  the 

400  Hz  servo  bandwidth.  An  amplitude  limitation  was  also  included  in  the 

simulation  to  keep  the  return  signal  from  exceeding  an  absolute  value  of  unity 

just  before  entering  the  loop  gain  amplifier  (noted  as  Gj  in  Fig.  37).  This 

modification,  which  avoids  instabilities  resulting  from  the  automatic  gain 

control,  is  used  in  the  laboratory  m del  of  the  COAT  array.  Simulations 

were  run  at  three  values  of  broad-band  signal-to-noise  ratio  (receiver  shot 

noise);  10^°,  40,  and  20.  When  (S/N)^^^^  = and  = 0 (no  speckle 

noise),  the  simulation  produces  a time  history  of  the  power  on  the  target 

glint  as  shown  in  Fig.  3 9(a).  When  the  speckle  noise  modulation  shown  in 

Fig.  39(c)  is  present  (a  = 0.  008);  the  time  history  in  Fig.  38(b)  results. 

rricLx 

In  this  case,  the  rms  value  of  is  only  0.  064  (around  = 0.  50),  but  the 
maximum  peak-to-peak  value  is  about  0.  4. 

The  behavior  shown  in  Fig.  39(b)  for  a _ ^ 0 is  typical.  In  general, 

the  speckle  noise  has  two  effects:  (1)  the  maximum  percent  convergence  level 

is  lowered,  ?nd  (2)  a random  oscillation  of  the  intensity  occurs  between  the 
new  maximum  convergence  level  and  a minimum  level.  This  oscillation  is 
much  slower  than  the  dither  oscillations  and  its  magnitude  varies  with  con- 
vergence level.  Figure  40  shows  a plot  of  the  maxima  and  minima  recorded 
during  these  simulations  as  a function  of  a . The  open  loop  gain,  G , was 
the  same  for  each  case.  At  high  average  convergence  levels  (low  values  of 

a ).  the  random  oscillations  have  much  smaller  magnitudes  than  at  the 
max  ’ 

low  convergence  levels  (high  values  of  a ). 

ms-x 

3 . Analysis  Compared  with  Computer  Simulation 

When  we  examine  other  speckle  power  curves  in  Fig.  36 

which  correspond  to  different  values  of  a and  cross  check  these  with  the 
^ max 

computer  simulation  data,  it  becomes  apparent  that  the  power  ratio  is  a very 
important  parameter,  defined  as 


Then  as  the  convergence  process  begins,  the  value  of  p is  much  greater  than 
1.  As  the  convergence  level  increases,  the  value  of  p gets  progressively 
smaller.  For  a given  a , the  COAT  system  should  converge  to  the  level 
corresponding  to  p = 1 (where  the  '^Pg>  and  curves  in  Fig.  35  cross). 

The  final  convergence  level  corresponding  to  p = I and  that  corresponding 
to  some  smaller  value  of  p. 

We  can  use  the  data  in  Fig.  40  along  with  the  analysis  developed 
above  to  verify  our  hypothesis  that  the  convergence  level  oscillates  between 
two  fixed  values  of  p and  to  find  these  values.  This  is  done  by  first  choosing 

p (e.  g.  , p = 1).  Equation  (12)  is  then  used  to  calculate  <Pp.>  as  a function 

2 1/2  ^ 

of  a (or  equivalently,  <I  > ).  Equations  (20)  and  (21)  are  then  used  to 

2 1/2 

calculate  a as  a function  of  <I  > . Using  this  method,  we  have  found 

mcLX  m ® 

that  0 =1.0  and  p = 0.  3 correspond  to  the  maximum  and  minimum 

^max  '^min 

convergence  levels.  Figure  41  compares  the  results  of  our  analytical 
calculations  with  the  computer -gene rated  data  of  Fig,  39.  The  agreement 
between  the  analytical  theory  and  the  simulation  data  is  excellent  for  signal- 
to-noise  (shot  noise)  ratios,  S/N,  between  10^^  and  40.  The  agreement  is 
not  as  good  for  the  S/n  = 20  case,  but  this  is  not  unexpected  since  the  com- 
puter simulations  have  always  shown  degradation  in  COAT  system,  perfor- 
mance when  S/N  was  chosen  less  than  40,  with  or  without  speckle  modulations. 
These  data  thus  appear  to  verify  the  original  hypothesis  that  motivated  much 
of  the  statistical  analysis,  namely  that  the  maximum  and  minimum  conver- 
gence levels  correspond  to  constant  values  of  the  dither-to- speckle  power 
ratio. 


The  good  agreement  between  our  simplified  analysis  and  the  computer 
simulation  data  has  given  us  confidence  that  we  can  accurately  predict  the 
performance  of  a multidither  COAT  system  in  the  presence  of  spurious 
modulations.  Our  work  to  date  has  been  with  a particular  COAT  system 
and  a hypothetical  model  for  the  spurious  modulations  (chosen,  however,  to 
present  maximum  difficulty  to  the  COAT  system).  We  expect  that  the  analy- 
sis can  be  extended  in  a straightforward  manner  so  that  it  is  independent  of 
the  speckle  modulation  model,  depending  only  on  the  statistics  of  the  modula- 
tion spectrum.  A contract  that  is  now  in  progress  (F30602 -76-C- 002 1 ) will 
accomplish  this  extension.  We  also  expect  that  modifications  to  the  servo. 
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Fig.  40,  COAT  convergence  levels  as  a function  of  speckle  noise  ampli- 
tude as  observed  by  computer  simulation. 
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Fig.  41.  Convergence  level  versus  amax  theory  versus  experiment 
for  an  18-element  COAT  array  (4^  = Z0°) 
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which  are  not  accounted  for  in  the  analysis  (the  details  of  the  low-pass  filter, 
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for  example),  will  affect  the  values  of  and  and  the  curves  <I^> 

versus  a for  a given  p.  The  ability  of  the  analysis  to  predict  the  con- 
max  ® 

vergence  levels  as  a function  of  a^^  (or  some  other  measure  of  the  speckle 
noise)  should  be  unaffected. 

At  this  point,  it  appears  that  the  speckle  modulation  problem  can  be 
viewed  as  a signal-to-noise  (S/N)  problem.  In  addition,  the  changes  in  COAT 
system  performance  caused  by  speckle  "noise"  are  independent  of  broadband. 
The  signal-to-noise  (shot,  thermal,  etc.  ) ratio  as  long  as  this  latter  S/N 
ratio  is  sufficient  for  stable  operation  the  absence  of  speckle  modulations. 

Tlie  COAT  system  signal-to- speckle  noise  ratio  varies  dramatically  with 
convergence  level,  however,  since  the  speckle  noise  within  the  ditherband 
competes  so  strongly  with  the  CO AT-generated  convergence  signals.  With 
this  viewpoint,  any  changes  in  the  COAT  servo  that  act  to  improve  the  S/N 
will  reduce  the  deleterious  effects  of  speckle  modulations.  For  example, 
we  expect  that  increasing  the  dither  amnlitude  or  the  amount  of  low-pass 
filtering  will  allow  operation  at  higher  speckle-noise  levels. 


V. 


SUMMARY  AND  CONCLUSIONS 


A.  Turbulence  and  Blooming  Compensation 

This  program  has  continued  to  build  the  experimental  capabilities 
of  the  DARPA/RADC  COAT  experimental  model  and  to  increase  its  ver- 
satility as  a breadb  >ird  for  studying  multidither  COAT  concepts  and  per- 
formance. This  system  has  proved  to  be  the  single  most  valuable  tool 
available  to  date  for  such  studies,  rivaled  only  by  the  detailed  computer 
simulations  that  have  been  developed  along  with  it.  The  COAT  system  now 
has  18  phase  control  channels,  two  tracking  control  channels,  and  one 
focus  control  channel.  For  propagation  studies,  there  are  several  thermal 
blooming  and  turbulence  generators,  including  a 100  m outdoor  path  at 
Hughef  Research  Laboratories. 

The  DARPA/RADC  COAT  system  has  demonstrated  that  an  improve- 
ment factor  of  about  1.  5 is  possible  against  blooming  distortions,  with  the 
compensation  occurring  only  for  distortions  in  the  first  30%  of  the  propaga- 
tion path.  Compensation  for  turbulence  and  fixed  distortions  is  not 
impaired  by  the  presence  of  thermal  blooming,  even  though  little  correction 
may  occur  for  the  blooming  itself.  This  observation  has  been  verified  by 
both  experiment  and  computer  simulation.  These  results,  along  with  those 
from  an  earlier  program,  indicate  that  the  achievement  of  the  high 
focal-plane  irradiance  should  not  be  limited  by  atmospheric  propagation, 
provided  that  thermal  blooming  distortions  are  minimized  and  that  a suit- 
able COAT  system  is  used. 

Computer  simulation  data  indicate  that  a return-wave  COAT  system 
may  not  be  a suitable  concept  when  thermal  blooming  is  present.  Such 
systems  can  remove  turbulence  and  fixed  distortions  just  as  a multidither 
system  does,  but  may  become  unstable  and  reduce  the  focal  plane 
irradiance  if  the  blooming  becomes  too  severe. 

In  all  of  the  experiments  performed  to  date,  the  multidither  tracking 
and  focus  (BAT)  controls  have  not  increased  the  average  peak  beam  irradi- 
ance significantly  over  that  produced  by  the  phase  controls  alone.  The 
only  exception  to  this  observation  occurred  with  artificial  turbulence 


distortions  introduced  with  n phase  screen.  With  a different  artificial 
phase  screen,  the  BAT  controls  had  no  effect.  This  somewhat  unexpected 
observation  may  be  a result  of  the  characteristics  of  our  artificial  tur- 
bulence screen  since  real  turbulence  should  have  significant  tilt  compo- 
nents (see  the  results  in  Section  IV  as  well  as  Ref.  11). 

B.  Complex  Target  Effects  Caused  by  Speckle  Modulations 

Problems  with  COAT  system  convergence  caused  by  speckle- 
modulations  have  not  been  observed  with  an  extended,  scotchlite  target 
used  with  the  18-channel  COAT  system.  Artificially-generated  amplitude 
modulations  have  affected  the  COAT  convergence  level,  however,  in  both 
experiments  and  computer  simulations.  It  remains  to  be  determined  whether 
real  targets  can  produce  significant  speckle  modulations  that  will  affect 
real  multidither  COAT  systems.  An  anal/tical  model  has  been  developed 
that  can  accurately  predict  the  convergence  level  and  stability  of  a multi- 
dither COAT  system  once  the  system  is  defined  and  once  the  statistics 
(such  as  the  power  spectral  density)  of  the  spurious  modulations  are  known. 
Additional  work  is  continuing  on  this  subject,  supported  by  DARPA/RADC 
Contract  F30602-76-C-0021. 

C.  Recommendations  for  Additional  Investigations 

The  work  on  this  contract  has  satisfied  the  original  goals  of  the 
program  and  has  contributed  substantially  to  the  understanding  of  COAT 
systems,  their  capabilities,  and  their  limitations.  There  are  several 
areas,  however,  where  additional  work  is  needed.  First,  we  need  to  deter- 
mine what  targets  and  scenarios  can  produce  speckle  modulations  that  are 
strong  enough  to  degrade  multidither  COAT  performance.  This  work  will 
involve  c.nalytical,  computer  simulation,  and  experimental  efforts  and 
should  inc.iude  some  experimental  work  at  10.6  pm,  or  another  convenient 
infrared  wavelength.  A current  contract,  F30602-76-C-00Z1,  ir  aimed 
at  some  of  this  work,  but  does  not  include  10.6  p.)ii  experimer  ral 
measurements . 
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A second  area  of  investigation  is  in  thermal  blooming  compensation. 
All  of  the  work  to  date,  at  Hughes  as  well  as  at  other  laboratories,  indi- 
cates that  present  COAT  techniques  will  not  do  much  for  thermal  blooming 
with  cw  lasers.  Other  techniques  such  as  different  multidither  concepts, 
irradiance  tailoring,  wavelength  selection,  and  pulsed  operation  should  be 
pursued  to  find  those  concepts  and  systems  that  minimize  thermal  bloomine 
or  increase  the  ability  of  COAT  to  compensate  for  those  blooming  dis- 
torHons  that  do  occur. 

There  is  still  much  to  be  done  in  the  area  of  computer  simulation. 

The  codes  developed  to  date  for  both  COAT  and  propagation  have  been  invalu- 
able in  identifying  and  studying  propagation  phenomena  and  compensation 
for  them.  They  have  also  been  essential  for  designing  and  evaluating  COAT 
hardware.  Additional  work  is  required  to  improve  these  codes  as  well  as 
to  utilize  them  in  studying  COAT  compensation  with  various  system  and 
target  scenarios.  We  need  to  develop  efficient  codes  that  can  simulate 
thermal  blooming  compensation  with  multidither  COAT  system  and  realistic, 
multiple-glint  targets. 

Since  most  of  the  high-energy  laser  applications  of  COAT  are  at 
infrared  wavelengths,  much  experimental  work  needs  to  be  done  to  demon- 
strate the  required  hardware  technology,  which  can  be  quite  different  with 
high-power,  infrared  lasers  than  with  low-power,  visible  wavelength  devices. 
The  HICLAS  program  (contract  N60921-76-C-0008  and  its  follow-on)  will 
address  some  of  these  technology  questions.  In  addition,  low-power  COAT 
experiments  with  infrared  wavelengths  are  in  order  to  investigate  target 
effects  such  as  speckle  and  glint  destruction  that  may  be  quite  different  at 
infrared  wavelengths  than  in  the  visible. 

Finally,  as  promising  new  COAT  servo  and/or  system  concepts  are 
invented,  they  should  be  investigated  in  low-power,  breadboard  experiments 
to  establish  what  advantages  they  may  have  in  compensation  performance, 
versatility,  simplicity  in  hardware  design  and  implementation,  or  relia- 
bility of  operation  under  a variety  of  operational  conditions. 
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MISSION 

of 

Rome  Air  Develojrment  Center 


HADC  is  the  principal  AFSC  organization  charged  with 
planning  and  executing  the  USAF  exploratory  and  advanced 
development  programs  for  information  sciences , intelli- 
gence, command,  control  and  communications  technology , 
products  and  services  oriented  to  the  needs  of  the  USAF . 
Primary  RADC  mission  areas  are  communications , electro- 
magnetic guidance  and  control,  surveillance  of  ground 
and  aerospace  objects,  intelligence  data  collection  and 
handling,  information  system  technology , and  electronic 
reliability,  maintainability  and  compatibility . RADC 
has  mission  responsibility  as  assigned  by  AFSC  for  de- 
monstration and  acquisition  of  selected  subsystems  and 
systems  in  the  intelligence , mapping,  charting , command, 
control  and  communications  areas. 


